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Preface 


View the promotional video on YouTube 


These are the lecture notes for my online Coursera course, Vector Calculus for Engineers. Students 
who take this course are expected to already know single-variable differential and integral calculus to 
the level of an introductory college calculus course. Students should also be familiar with matrices, 
and be able to compute a three-by-three determinant. 

I have divided these notes into chapters called Lectures, with each Lecture corresponding to a 
video on Coursera. I have also uploaded all my Coursera videos to YouTube, and links are placed at 
the top of each Lecture. 

There are some problems at the end of each lecture chapter. These problems are designed to 
exemplify the main ideas of the lecture. Students taking a formal university course in multivariable 
calculus will usually be assigned many more problems, some of them quite difficult, but here I follow 
the philosophy that less is more. I give enough problems for students to solidify their understanding 
of the material, but not so many that students feel overwhelmed. I do encourage students to attempt 
the given problems, but, if they get stuck, full solutions can be found in the Appendix. I have also 
included practice quizzes as an additional source of problems, with solutions also given. 


JEFFREY R. CHASNOV 
Hong Kong 
October 2019 
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Week I 


Vectors 


In this week's lectures, we learn about vectors. Vectors are line segments with both length and direction, and 
are fundamental to engineering mathematics. We will define vectors, how to add and subtract them, and how 
to multiply them using the scalar and vector products (dot and cross products). We use vectors to learn some 
analytical geometry of lines and planes, and introduce the Kronecker delta and the Levi-Civita symbol to prove 
vector identities. The important concepts of scalar and vector fields are discussed. 


Lecture 1 | Vectors 


View this lecture on YouTube 


We define a vector in three-dimensional Euclidean space as having a length (or magnitude) and a 
direction. A vector is depicted as an arrow starting at one point in space and ending at another point. 
All vectors that have the same length and point in the same direction are considered equal, no matter 
where they are located in space. (Variables that are vectors will be denoted in print by boldface, and 
in hand by an arrow drawn over the symbol.) In contrast, scalars have magnitude but no direction. 
Zero can either be a scalar or a vector and has zero magnitude. The negative of a vector reverses its 
direction. Examples of vectors are velocity and acceleration; examples of scalars are mass and charge. 

Vectors can be added to each other and multiplied by scalars. A simple example is a mass m acted 
on by two forces F} and Fy. Newton's equation then takes the form ma = Fi + Fy, where a is the 


acceleration vector of the mass. Vector addition is commutative and associative: 
A+B=B+A, (A+B)+C=A+(B+C); 
and scalar multiplication is distributive: 
k(A+B)=kA+kB. 


Multiplication of a vector by a positive scalar changes the length of the vector but not its direction. 
Vector addition can be represented graphically by placing the tail of one of the vectors on the head of 
the other. Vector subtraction adds the first vector to the negative of the second. Notice that when the 
tail of A and B are placed at the same point, the vector B — A points from the head of A to the head 
of B, or equivalently, the tail of — A to the head of B. 


B- A 
A+B=-B+A A 

B 
Vector addition Vector subtraction 


4 LECTURE 1. VECTORS 


Problems for Lecture 1 


1. Show graphically that vector addition is associative, that is, (A+ B) + C = A + (B + C). 


2. Using vectors, prove that the line segment joining the midpoints of two sides of a triangle is parallel 
to the third side and half its length. 


Solutions to the Problems 


Lecture 2 Cartesian coordinates 


View this lecture on YouTube 


To solve a physical problem, we usually impose a coordinate system. The familiar three-dimensional 
x-y-z coordinate system is called the Cartesian coordinate system. Three mutually perpendicular 
lines called axes intersect at a point called the origin, denoted as (0,0,0). All other points in three- 
dimensional space are identified by their coordinates as (x,y,z) in the standard way. The positive 
directions of the axes are usually chosen to form a right-handed coordinate system. When one points 
the right hand in the direction of the positive x-axis, and curls the fingers in the direction of the 
positive y-axis, the thumb should point in the direction of the positive z-axis. 

A vector has a length and a direction. If we impose a Cartesian coordinate system and place the 
tail of a vector at the origin, then the head points to a specific point. For example, if the vector A has 
head pointing to (A1, A2, A3), we say that the x-component of A is Ay, the y-component is A», and 
the z-component is A5. The length of the vector A, denoted by |A|, is a scalar and is independent of 
the orientation of the coordinate system. Application of the Pythagorean theorem in three dimensions 


|A| = y AZ +43 + A}. 


We can define standard unit vectors i, j and k, to be vectors of length one that point along the 


results in 


positive directions of the x-, y-, and z-coordinate axes, respectively. Using these unit vectors, we write 
A = Aıi+ Aoj + Ask. 


With also B = Bi + B2j + B3k, vector addition and scalar multiplication can be expressed component- 
wise and is given by 


A+ B = (A1 +B1) + (A5 + B5)j + (Aa + B3), cA=cAıi+cA2j +cA3k. 


The position vector, r, is defined as the vector that points from the origin to the point (x, y, z), and is 


used to locate a specific point in space. It can be written in terms of the standard unit vectors as 
T — Xi -- yj t zk. 


A displacement vector is the difference between two position vectors. For position vectors rı and rp, 
the displacement vector that points from the head of rı to the head of rz is given by 


r2 — 11 = (x2 — x1) + (y — y1)j + 22 — zi). 


LECTURE 2. CARTESIAN COORDINATES 


Problems for Lecture 2 


a) Given a Cartesian coordinate system with standard unit vectors i, j, and k, let the mass mı be 
at position r1 = x11 + y1j + zik and the mass m» be at position ra = x» + y2j + zak. In terms 
of the standard unit vectors, determine the unit vector that points from mı to mn. 


Newton's law of universal gravitation states that two point masses attract each other along the 
line connecting them, with a force proportional to the product of their masses and inversely 
proportional to the square of the distance between them. The magnitude of the force acting on 


each mass is therefore 
mm; 


F=G E 
r2 


where mı and m» are the two masses, r is the distance between them, and G is the gravitational 
constant. Let the masses m] and mz be located at the position vectors rı and r2. Write down the 
vector form for the force acting on mı due to its gravitational attraction to mp. 


Solutions to the Problems 


Lecture 3 | Dot product 


View this lecture on YouTube 


We define the dot product (or scalar product) between two vectors A = Ai + A2j + A3k and 
B = Byi+ Boj + B3k as 
A- B= A,B, + A2B2 + A3B3. 


One can prove that the dot product is commutative, distributive over addition, and associative with 


respect to scalar multiplication; that is, 
A:-B=B-A, A:-(B+C)=A:-B+A:C, A-(cB)=(cA)-B=c(A-B). 


A geometric interpretation of the dot product is also possible. Given any two vectors A and B, 
place the vectors tail-to-tail, and impose a coordinate system with origin at the tails such that A is 
parallel to the x-axis and B lies in the x-y plane, as shown in the figure. The angle between the two 
vectors is denoted as 0. 

Then in this coordinate 
system, A = |Ali, B = 
|B| cos + | B| sin 05, and 


A-B=|A|B|cos$, 


a result independent of the 
choice of coordinate system. 
If A and B are parallel, then B 
0 =0and A- B = |A||B| 
and in particular, A- A = 
|A|?. If A and B are per- 
pendicular, then 9 = 7/2 
and A. B =0. 0 


|B| cos 0 A 


8 LECTURE 3. DOT PRODUCT 


Problems for Lecture 3 
1. Using the definition of the dot product A - B = A,B, + A2B2 + A3B3, prove that 
a) A-B=B-A; 
b) A-(B+C)=A-B+A-C; 
c) A- (kB) = (kA): B =k(A- B). 
2. Determine all the combinations of dot products between the standard unit vectors i, j, and k. 
3. Let C = A — B. Calculate the dot product of C with itself and thus derive the law of cosines. 


Solutions to the Problems 


Lecture 4 | Cross product 


View this lecture on YouTube 


We define the cross product (or vector product) between two vectors A = Ai + A53 + Ask and 
B = Bi + Boj + B3k as 


i o4 E 
AxB=|Aı Az A3| = (A2B3 — A3B2)i + (43B1 — A1B3)j + (A1Bz — A2B1)k. 
By By Bs 


Defining the cross product by a determinant serves to remember the formula. One can prove that 
the cross product is anticommutative, distributive over addition, and associative with respect to scalar 


multiplication; that is 


Ax B—-BxA, Ax(B+C)=AxB+AxC, A x (cB) = (cA) x B=c(Ax B). 


A axb 


B |B|sin 0 S £ 


a Right-hand rule (from Wikipedia) 


A geometric interpretation of the cross product is also possible. Given two vectors A and B with 
angle 0 between them, impose a coordinate system so that A is parallel to the x-axis and B lies in the 
x-y plane. Then A = |Ali, B = |B|cos6i + | B| sin0j, and A x B = |A||B|sin@k. The coordinate- 
independent relationship is 

|A x B| = |A||B| sin 6, 


where 0 lies between zero and 180°. Furthermore, the vector A x B points in the direction perpendic- 
ular to the plane formed by A and B, and its sign is determined by the right-hand rule. Also, observe 
that |A x B| is the area of the parallelogram whose adjacent sides are the vectors A and B. 


10 LECTURE 4. CROSS PRODUCT 


Problems for Lecture 4 
1. Using properties of the determinant, prove that 
a Ax B=-B x A; 
b) Ax (B-C)- Ax B- AxC; 
c) A x (kB) = (kA) x B —k(A x B). 
2. Determine all the combinations of cross products between the standard unit vectors i, j, and k. 


3. Show that the cross product is not in general associative. That is, find an example using unit vectors 
such that 
Ax(BxC)z(Ax B)xC. 


Solutions to the Problems 


Practice Quiz | Vectors 


1. Let A, B and C be any vectors. Which of the following statements is false? 
a A-B=B-A 
b) A+(B+C)=(A+B)+C 
c) Ax(BxC)=(AxB)xC 
d A-(B+C)=A-B+A-C 
2. Let A = ai + mmj +. azk and B = b1i +b2j +b3k. Then (A x B)- j is equal to 
a) a2b3 — aab; 
b) azb¡ — a1b3 
c) aybz — a5bi 
d) a1b3 — azbı 
3. Which vector is not equal to zero? 
a) ix (j x k) 
b) (ixj)xk 
c) (ixi)xj 
d) à x (ixj) 


Solutions to the Practice quiz 
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PRACTICE QUIZ. VECTORS 


Lecture5 | Analytic geometry of lines 


View this lecture on YouTube 


In two dimensions, the equation for a line in slope-intercept form is y = mx + b, and in point-slope 
form is y — y; — m(x — x1). In three dimensions, a line is most commonly expressed as a parametric 
equation. 

Suppose that a line passes through a point with position vector ro and in a direction parallel to 
the vector u. Then, from the definition of vector addition, we can specify the position vector r for any 
point on the line by 

T — TQ 4 ut, 


where t is a parameter that can take on any real value. 

This parametric equation for a line has clear physical meaning. If r is the position vector of a 
particle, then u is the velocity vector, and t is the time. In particular, differentiating r — r(t) with 
respect to time results in dr/dt = u. 

A nonparametric equation for the line can be obtained by eliminating t from the equations for the 


components. The component equations are 
X = xo + ut, y = yo + uot, Z = Zo + ust; 


and eliminating t results in 


Xx— X0 _Y— Y0 _2—20 
u uz uz ` 


Example: Find the parametric equation for a line that passes through the points (1,2,3) and (3,2,1). Determine 
the intersection point of the line with the z = 0 plane. 


To find a vector parallel to the direction of the line, we first compute the displacement vector between 


the two given points: 


u = (3—1)i + (2—2)j + (1—3)k = 2i — 2k. 


Choosing a point on the line with position vector rg = i + 2j + 3k, the parametric equation for the 
line is given by 


r = ro + ut =i+2j+3k+t(2i— 2k) = (1 +2t)i+2j + (3 — 2t)k. 


The line crosses the z = 0 plane when 3 — 2t = 0, or t = 3/2, and (x,y) = (4,2). 
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14 LECTURE 5. ANALYTIC GEOMETRY OF LINES 


Problems for Lecture 5 


1. Find the parametric equation for a line that passes through the points (1,1,1) and (2,3,2). Deter- 


mine the intersection point of the line with the x — 0 plane, y — 0 plane and z — 0 plane. 


Solutions to the Problems 


Lecture 6 | Analytic geometry of planes 


View this lecture on YouTube 


A plane in three-dimensional space is determined by three non-collinear points. Two linearly indepen- 
dent displacement vectors with direction parallel to the plane can be formed from these three points, 
and the cross-product of these two displacement vectors will be a vector that is orthogonal to the 
plane. We can use the dot product to express this orthogonality. 

So let three points that define a plane be located by the position vectors r1, r2, and r3, and construct 
any two displacement vectors, such as s, = r2 — rı and s? = r3 — r2. A vector normal to the plane is 
given by N = s, x 82, and for any point in the plane with position vector r, and for any one of the 
given position vectors r;, we have N - (r — rj) = 0. With r =xi+tyj+zk, N = ai 4 bj +ck and 
d = N - r;, the equation for the plane can be written as N -r = N «rj, or 


ax + by -- cz = d. 
Notice that the coefficients of x, y and z are the components of the normal vector to the plane. 


Example: Find an equation for the plane defined by the three points (2,1,1), (1,2,1), and (1,1,2). Determine 
the equation for the line in the x-y plane formed by the intersection of this plane with the z = 0 plane. 


To find two vectors parallel to the plane, we compute two displacement vectors from the three points: 


s=(1-2Di+(2-Dji+(1-Dk=-4+4  s-—(-1ic-(1-2)j + 2—1)k- —j 4k. 


A normal vector to the plane is then found from 
i k 
N = sı X s2 = |—1 1 0|=:+Jj+k. 
0 —1 1 
And the equation for the plane can be found from N -r = N - r4, or 


(i+j+k) (xityj+zk)=(i+j+k) (i+j+k), or x+y+z=4. 


The intersection of this plane with the z = 0 plane forms the line given by y — —x +4. 
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16 LECTURE 6. ANALYTIC GEOMETRY OF PLANES 


Problems for Lecture 6 


1. Find an equation for the plane defined by the three points (—1, —1, —1), (1, 1,1) and (1, —1,0). 
Determine the equation for the line in the x-y plane formed by the intersection of this plane with the 


z — 0 plane. 


Solutions to the Problems 


Practice Quiz | Analytic geometry 


1. The line that passes through the points (0,1,1) and (1,0, —1) has parametric equation given by 


a) ti+(1+0j+(1420)k 


b) ti+ (1— O3 + (13-2£)k 


c) ti - (14-£)j  (1—20k 


d) i - (1—t)j + (1 -20)k 


2. The line of Question 1 intersects the z — 0 plane at the point 


3. The equation for the line in the x-y plane formed by the intersection of the plane defined by the 
points (1,1, 1), (1,1,2) and (2, 1,1) and the z = 0 plane is given by 


a)y=x 
b)y=x+1 
c)y=x-1 
d) y —1 


Solutions to the Practice quiz 
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PRACTICE QUIZ. ANALYTIC GEOMETRY 


Lecture 7 | Kronecker delta and 
Levi-Civita symbol 


View this lecture on YouTube 


We define the Kronecker delta bij to be +1 if i = j and 0 otherwise, and the Levi-Civita symbol 
Eijk to be +1 if i, j, and k are a cyclic permutation of (1,2,3) (that is, one of (1,2,3), (2,3,1) or (3,1,2)); 
—1 if an anticyclic permutation of (1,2,3) (that is, one of (3,2,1), (2,1,3) or (1,3,2)); and 0 if any two 


indices are equal. More formally, 
1, ifi=j; 
ij = o 
0, ifizj, 


+1, if (i,j,k) is (1,2,3), (2,3,1) or (3,1,2); 
ex = 4 —1, if (i,j,k) is (3,2,1), (2,1,3) or (1,3,2); 
0, ifi=jorj=kork=1. 


For convenience, we will use the Einstein summation convention when working with these sym- 
bols, where a repeated index implies summation over that index. For example, we will write an 
expression such as Pd A;B; more simply as A;B;. In a product of components, any index can be 
repeated at most one time, though multiple indices may be repeated. 

Some interesting examples using the Kronecker delta and the Levi-Civita symbols are 


"m =3 and Cijk€ijk = 6, 


where ö;; = 911 + 022 + 633 and €ijk€ijk implies a sum over i, j, and k, containing a total of 8? — 27 terms 
in the sum, where six of the terms are equal to one and the remaining terms are equal to zero. 

We note here a remarkable relationship between the product of Levi-Civita symbols and the Kro- 
necker delta, given by the determinant 


Ôi Sim Sin 
Cijk€lmn = dil dim din = Oi (SjmOkn x ÓjinÓkm) — Sim (jk — 951041) + Óin (Sjrökm — Oim0%1)- 
Ód Okm kn 
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20 LECTURE 7. KRONECKER DELTA AND LEVI-CIVITA SYMBOL 


Problems for Lecture 7 
1. Prove the following cyclic and anticyclic permutation identities: 
A) Eijk = €jki = €ki; 
b) Eijk = —€jikr Eijk = —Ekjir Eijk = —Eikj- 
2. Verify the cross-product relation [A x B]; = €jjxA¡Bx by considering i = 1,2,3. 
3. Prove the following Kronecker-delta identities: 
a) dj Aj = Ai; 
b) óióyj = bij. 
4. Given the most general identity relating the Levi-Civita symbol to the Kronecker delta, 
€ijk€Imn = Óil(ÓjmÓkn — Ojndkm) — Fim(Oj15kn — Ójn9x1) + Sin (Oj Ókm — Ójm0r1 ), 
prove the following simpler and more useful relations: 
A) €iik€imn = ÓjmOkn — OjnÓkm; 
b) éijk€ijn = 2ôkn- 


Solutions to the Problems 


Lecture 8 | Vector identities 


View this lecture on YouTube 


Four sometimes useful vector identities are 


(scalar triple product) A:(BxC)- B.(Cx A) 2 C. (Ax B), 

(vector triple product) Ax(BxC)=(A-C)B-(A-B)C 

(scalar quadruple product) (Ax B)-(CxD)=(A-C)(B-D)-(A-D)(B-C) 
(vector quadruple product) (Ax B)x(CxD)=((AxB)-D)C-((Ax B): C)D 


Parentheses are optional when expressions have only one possible interpretation, but for clarity they 
are often written. The first and third identities equate scalars whereas the second and fourth identities 
equate vectors. To prove two vectors A and B are equal, one must prove that their components are 
equal, that is A; = Bj for an arbitrary index i. 

Proofs of these vector identities can make use of the Kronecker delta, the Levi-Civita symbol and 
the Einstein summation convention, where repeated indices are summed over. The algebraic toolbox 
that we will need contains the cyclic permutation identities for the Levi-Civita symbol: 


Eijk = €jki = €kijo 
the contraction of two Levi-Civita symbols: 
Eijk€imn = ÓjmÓkn — ÓjnÓkm; 
the contraction of the Kronecker delta with a vector: 
dijA j= Aj; 
and the scalar and vector products written as 


A-B= A;B;, [A x B]; — Ejjk Aj By. 


We will use the notation [...]; to mean the ith component of the bracketed vector. 
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22 LECTURE 8. VECTOR IDENTITIES 


Problems for Lecture 8 


1. Remove all optional parentheses from the four vector identities 
(scalar triple product) A:(BxC)-—- B.(Cx A) 2 C. (Ax B), 
(vector triple product) Ax(BxC)=(A-C)B-(A-B)C 
(scalar quadruple product) (Ax B)-(C x D)=(A-C)(B-D)-(A:D)(B:C) 
(vector quadruple product) (Ax B) x (Cx D) 2 ((Ax B): D)C—((Ax B): C)D 


Solutions to the Problems 


Lecture 9 | Scalar triple product 


View this lecture on YouTube 


The scalar triple product, which can be written as either A - (B x C) or A: B x C, satisfies 
A-(BxC)— B-.(Cx A) - C (Ax B), 
that is, its value is unchanged under a cyclic permutation of the three vectors. This identity can be 


proved using the cyclic property of the Levi-Civita symbol: 
A: (B x C) = Ai€ijkBjCk = BjeikiCkAi =B. (C x A), 


and 
B. (C x A) = Bje ki Ck Ai = CreriAiBj = C> (A x B). 


We can also write the scalar triple product as a three-by-three determinant. We have 
j k A, Az Ag 

Bo B3 = B4 B5 B3 . 

Cy & C3 


i 
Ci GQ C3 


The absolute value of a three-by-three determinant, and therefore that of the scalar triple product, is 
the volume of the parallelepiped defined by the three row vectors A, B and C, as shown in the figure 


below. 
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24 LECTURE 9. SCALAR TRIPLE PRODUCT 


Problems for Lecture 9 


1. Consider the scalar triple product A - (B x C). Prove that if any two vectors are equal, then the 


scalar triple product is zero. 


2. Show that swapping the positions of the operators without re-ordering the vectors leaves the scalar 
triple product unchanged, that is, 
A-BxC=AxB-C. 


3. It is sometimes useful to define a notation where the unit vectors are distinguished by their index. 


That is, we define e4 = i, e? = j and e3 = k. Prove that 
ei (ej X ex) = Eijk 


where ejj, is the usual Levi-Civita symbol. 


Solutions to the Problems 


Lecture 10 | Vector triple product 


View this lecture on YouTube 


The vector triple product satisfies 
Ax(BxC)-—(A-C)B -(A- B)C. 


This identity can be proved using Kronecker delta and Levi-Civita symbol identities as follows: 


[A x (B x C)]i = eig Aj|B x C]x ([A x X]; = €jjx AjXx) 
= eic AjEktmBICm([B x C]k = €k1mBiCm) 
= €kij€kIm AjBICm (Eijk = €xij) 
= (Ói10jm — Sim j1)AjBiCm (€xij€kIm = Sijm — Sim9j1) 
= AjB|C; — A;Bj;Cjo; Bi = Bi, etc. 
= (A: C)B - (A- B)C];. (AjC; = A: C, etc.) 
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26 LECTURE 10. VECTOR TRIPLE PRODUCT 


Problems for Lecture 10 


1. Prove the Jacobi identity: 
Ax(BxC)-Bx(Cx A) -Cx(Ax B)=0. 
2. Prove that the scalar quadruple product satisfies 
(Ax B): (Cx D) (A. CB. D) - (A. D(B.C). 
3. Prove Lagrange's identity in three dimensions: 
¡Ax BP? =|AP|B?-(4- B. 

4. Prove using vector triple products that the vector quadruple product satisfies 

(Ax B) x (Cx D) = ((Ax B): D)C — ((Ax B): C)D. 


Solutions to the Problems 


Practice Quiz | Vector algebra 


1. The expression ej;;€j;,, is equal to 
A) dimOkn — OinOkm 
b) Sinókm — ÓjmÓkn 
C) Skmdit — Óx1Óim 
d) dk0im — Ökmõil 
2. The expression A x (B x C) is always equal to 
a Bx (Cx A) 
b) Ax (C x B) 
c) (Ax B)x C 
d) (Cx B)x A 
3. Which of the following expressions may not be zero? 
a) A- (B x B) 
b) A: (Ax B) 
c) Ax (A x B) 
d) B-(A x B) 


Solutions to the Practice quiz 
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PRACTICE QUIZ. VECTOR ALGEBRA 


Lecture 11 


View this lecture on YouTube 


In physics, scalars and vectors can be functions of both space and time. We call these types of functions 


fields. For example, the temperature in some region of space is a scalar field, and we can write 


where we use the common notation of a semicolon on the right-hand side to separate the space and 
time dependence. Notice that the position vector r is used to locate the temperature in space. As 


Scalar and vector fields 


T(r,t) = T(x, y zt), 


another example, the velocity vector u of a flowing fluid is a vector field, and we can write 


u(r,t) = um (x,y, z;t)i + uo(x,y,z; t)j +u3(x, y,2;t)k. 


The equations governing a field are called the field equations, and these equations commonly take 
the form of partial differential equations. For example, the equations for the electric and magnetic 
vector fields are the famous Maxwell's equations, and the equation for the fluid velocity vector field 
is called the Navier-Stokes equation. The equation for the scalar field (called the wave function) in 


non-relativistic quantum mechanics is called the Schródinger equation. 


There are many ways to visu- 
alize scalar and vector fields, and 
one tries to make plots as infor- 
mative as possible. On the right is 
a simple visualization of the two- 
dimensional vector field given by 


—yictxj 
Boxy) = We 


, 


where the vectors at each point 


are represented by arrows. 
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30 LECTURE 11. SCALAR AND VECTOR FIELDS 


Problems for Lecture 11 


1. List some physical examples of scalar and vector fields. 


Solutions to the Problems 


Week II 


Differentiation 


In this week's lectures, we learn about the derivatives of scalar and vector fields. We define the partial derivative 
and derive the method of least squares as a minimization problem. We learn how to use the chain rule for 
a function of several variables, and derive the triple product rule used in chemistry. From the del differential 
operator, we define the gradient, divergence, curl and Laplacian. We learn some useful vector derivative identities 
and how to derive them using the Kronecker delta and Levi-Civita symbol. Vector identities are then used to 
derive the electromagnetic wave equation from Maxwell’s equations in free space. Electromagnetic waves are 
fundamental to all modern communication technologies. 
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Lecture 12 | Partial derivatives 


View this lecture on YouTube 


For a function f — f(x,y) of two variables, we define the partial derivative of f with respect to x 


Of _ im thv) Sy) 
Ox  h>30 h ’ 


as 


and similarly for the partial derivative of f with respect to y. To take a partial derivative with respect 
to a variable, take the derivative with respect to that variable treating all other variables as constants. 
As an example, consider 

f(x,y) = 2:2? +y. 
We have 


Of pag Of a 2 
à: xy, y I 


Second derivatives are defined as the derivatives of the first derivatives, so we have 


Pf 2 Bf 3 
3x —12xy, ay = 4x” + 6y; 


and for continuous differentiable functions, the mixed second partial derivatives are independent of 


the order in which the derivatives are taken, 


of 2, f 
dxdy Hy dydx' 


To simplify notation, we introduce the standard subscript notation for partial derivatives, 


of of _ of 


of of 
fx = da fy = , Fix = 3x2’ fxy = axdy’ fuy dy?” 


ay etc. 


The Taylor series of f(x,y) about the origin is developed by expanding the function in a multivariable 
power series that agrees with the function value and all its partial derivatives at (x,y) = (0,0). We 


have 
1 
f(x,y) = f + fxx + fyy + 2! p + 2fxyXy + fu?) +.... 


The function and all its partial derivatives on the right-hand side are evaluated at (0,0) and are 
constants. By taking partial derivatives it is evident that f(0,0) = f, fx(0,0) = fx, fy(0,0) = fy, and 


so on as the infinite series continues. 
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34 LECTURE 12. PARTIAL DERIVATIVES 


Problems for Lecture 12 


1. Compute the three partial derivatives of 


1 
x2 ES y? + z2)” i 


fige) os 
2. Given the function f = f (t, x), find the Taylor series expansion of the expression 
f (t - & At, x + BAtf(t,x)) 


to first order in At. 


Solutions to the Problems 


Lecture 13 | The method of least squares 


View this lecture on YouTube 


Local maxima and minima of a multivariable function can be found by computing the zeros of the 

partial derivatives. These zeros are called the critical points of the function. A critical point need not 

be a maximum or minimum, for example it might be a minimum in one direction and a maximum in 

another (called a saddle point), but in many problems maxima or minima may be assumed to exist. 

Here, we illustrate the procedure for minimizing a function by solving the least-squares problem. 
Suppose there is some experimental data that you want 

to fit by a straight line (illustrated on the right). In general, 

let the data consist of a set of n points given by (x1,y1), ..., 


(Xn, Yn). Here, we assume that the x values are exact, and the 
y values are noisy. We further assume that the best fit line 
to the data takes the form y = fo + Bıx. Although we know > 
that the line can not go through all the data points, we can 


try to find the line that minimizes the sum of the squares of 


the vertical distances between the line and the points. 
Define this function of the sum of the squares to be 


n 


f (Bor b1) = } (Bo + bixi - i 


i=1 


Here, the data is assumed given and the unknowns are the fitting parameters ßo and f4. It should be 
clear from the problem specification, that there must be values of By and f that minimize the function 
f = f(Bo, B1). To determine, these values, we set df /0Bo = 9f /061 = 0. This results in the equations 


Y: (Bo + Bux — yi) =0, Yx (Bo + B1x; — yi) =0. 
i=1 


i=1 


We can write these equations as a linear system for Do and D, as 


n n n n n 
Bon+Bidixi=Yoyi, Bod mit Bi ox? = Ve xvi 
1=1 i=1 i=1 i=1 i=] 


The solution for Bo and B, in terms of the data is given by 


Bx XX 
nu — (Lx)? 


| nye xii - (xi) (Lyi) 
np -— (Lx)? 


Bo , Bi 


, 


where the summations are from i — 1 to n. 
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36 LECTURE 13. THE METHOD OF LEAST SQUARES 


Problems for Lecture 13 


1. By minimizing the sum of the squares of the vertical distance between the line and the points, 
determine the least-squares line through the data points (1,1), (2,3) and (3,2). 


Solutions to the Problems 


Lecture 14 | Chain rule 


View this lecture on YouTube 


Partial derivatives are used in applying the chain rule to a function of several variables. Consider 
a two-dimensional scalar field f = f(x,y), and define the total differential of f to be 


df = f(x t dx,y + dy) — f(x,y). 


We can write df as 


df = [f(x +dx,y + dy) — f(x,y + dy)] + [f(xy + dy) — f(x,y)] 


-fy f 
= 3, 7* + ay 


If one has f = f(x(t),y(t)), say, then division of df by dt results in 


df  Ofdx , Of dy 
dt  Oxdt Oydt 


And if one has f = f(x(r,0), y(r,0)), say, then the corresponding chain rule is given by 


of ofóx | of ay af ofóx | af oy 
or  Oxor Oyor' 90 0x90 | Oy 00 


dx 


d 
L u(t,x(t)). Determine a formula for E in terms of u and its 


Example: Consider the differential equation 


partial derivatives. 


Applying the chain rule, we have at time t, 


The above formula is called the material derivative and in three dimensions forms a part of the Navier- 


Stokes equation for fluid flow. 
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38 LECTURE 14. CHAIN RULE 


Problems for Lecture 14 


1. Let f(x, y) = e*", with x = rcos0 and y = rsin@. Compute the partial derivatives 9f /dr and df /00 
in two ways: 


a) Use the chain rule on f = f (x(r,0),y(r,0)); 
b) Eliminate x and y in favor of r and 0 and compute the partial derivatives directly. 


Solutions to the Problems 


Lecture 15 | Triple product rule 


View this lecture on YouTube 


Suppose that three variables x, y and z are related by the equation f(x,y,z) = 0, and that it is possible 
to write x = x(y,z), y = y(x,z) and z = z(x,y). Taking differentials of x and z, we have 


ox ox Oz Oz 
dx = z” + 3; 7^ dz — 3: 7* + a; 


We can make use of the second equation to eliminate dz in the first equation to obtain 


Ox ox / Əz Oz 
dx dy t Jz (Zax T E 


Ox oz ox | Oxoz 
(: ese) dy (H+ ES) ay, 


Since dx and dy can be independent variations, the terms in parentheses must be zero. The left-hand- 


or collecting terms, 


side results in the reciprocity relation 
Ox Oz _ 
azox ^ 
which states the intuitive result that 0z/dx and 0x/0z are multiplicative inverses of each other. The 


right-hand-side results in 
ox Ox Oz 


y ~~ Bz By’ 
which when making use of the reciprocity relation, yields the counterintuitive triple product rule, 


dx dy 02 


dy oz ox 
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LECTURE 15. TRIPLE PRODUCT RULE 


Lecture 16 | Triple product rule 
(example) 


View this lecture on YouTube 


Example: Demonstrate the triple product rule using the ideal gas law. 


The ideal gas law states that 
PV = nRT, 


where P is the pressure, V is the volume, T is the absolute temperature, n is the number of moles of 
the gas, and R is the ideal gas constant. We say P, V and T are the state variables, and the ideal gas 
law is a relation of the form 

f(P,V,T) = PV —nRT =0. 


We can write P = P(V,T), V = V(P, T) and T = T(P, V), that is, 


o nRT nRT PV. 


p V” P nR' 


and the partial derivatives are given by 


The triple product results in 
oP OVOT _ nRT nR VN nRT _ 1 
ƏV OT oP v2 P nR PV i 


where we make use of the ideal gas law in the last equality. 
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42 LECTURE 16. TRIPLE PRODUCT RULE (EXAMPLE) 


Problems for Lecture 16 


1. Suppose the three variables x, y and z are related by the linear expression ax + by + cz = 0. Show 
that x, y and z satisfy the triple product rule. 


2. Suppose the four variables x, y, z and t are related by the linear expression ax + by + cz + dt = 0. 
Determine a corresponding quadruple product rule for these variables. 


Solutions to the Problems 


Practice Quiz | Partial derivatives 


1 9 
1. Let f(x,y,z) = een The mixed second partial derivative xd is equal to 


2(x +y) 
x2 + y2 + 22)572 


(x+y)? 
x2 + y2 + z2)5/2 


42 
x2 + y2 + 22)5/2 


3xy 


a) x2 + y2 + 22)5/2 


2. The least-squares line through the data points (0,1), (1,3), (2,3) and (3,4) is given by 


7 9x 
OU AE 
5 9x 
M NT 
7 10x 
ui ge 
5 10x 
UT e ur 
3. Let f = f(x, y) with x = rcos0 and y = rsin 0. Which of the following is true? 
Əf _ WO i 
du x RT 
Əf A 
Jo cct 
OF _ 3f f 
Da gg tan 


of _ of, of 
a en 


Solutions to the Practice quiz 
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PRACTICE QUIZ. PARTIAL DERIVATIVES 


Lecture 17 | Gradient 


View this lecture on YouTube 


Consider the three-dimensional scalar field f — f(x, y, z), and the differential df, given by 


A 1, Sg 
df 3:4* | 2/9 H 3 12- 


Using the dot product, we can write this in vector form as 


ie (e | Ji | a) (Id EaR] YE ety 


where in Cartesian coordinates, 
_Əf. PET I 
Vf Jx TH ay JT E k 


is called the gradient of f. The nabla symbol V is pronounced "del" and V f is pronounced "del-f". 
Another useful way to view the gradient is to consider V as a vector differential operator which has 


the form 
0 d d 


= ET | E 


Because of the properties of the dot product, the differential df is maximum when the infinitesimal 
displacement vector dr is along the direction of the gradient V f. We then say that V f points in the 
direction of maximally increasing f, and the magnitude of V f gives the slope (or gradient) of f in that 
direction. 

Example: Compute the gradient of f (x,y,z) — xyz. 


The partial derivatives are easily calculated, and we have 


Vf=yzi+txzj+xyk. 
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46 LECTURE 17. GRADIENT 


Problems for Lecture 17 


1. Let the position vector be given by r = xi + yj + zk, with r = |r| = yx? +y? 4- zZ. Compute the 
gradient of the following scalar fields and write the result in terms of r and r. 
a) p(x, y, z) = x? +y +27; 
b) (uy) = = 
G S +y +z? 


Solutions to the Problems 


Lecture 18 | Divergence 


View this lecture on YouTube 


Consider in Cartesian coordinates the three-dimensional vector field, u = u (x,y,z)i + uo(x,y,z)3 + 
u3(x,y,z)k. The divergence of u, denoted as V - u and pronounced “del-dot-u”, is defined as the 
scalar field given by 


3,2 9 TM 
Vu e: Hi, | ka) - (ut + u2j +u3k) 


| ðu, , du, 
dx | dy | oO. 


Here, the dot product is used between a vector differential operator V and a vector field u. The diver- 
gence measures how much a vector field spreads out, or diverges, from a point. A more math-based 
description will be given later. 


Example: Let the position vector be given by r = xi + yj + zk. Find V r. 


A direct calculation gives 


ð ð 9 
V.r Jx dy? j 3:4 3. 
Example: Let F = „per all r #0. Find V - F. 
T 


Writing out the components of F, we have 


x y zZ 


Narren (x2 + yim (x2 + y2 + 223727 Pr rere” 


Using the quotient rule for the derivative, we have 


Or (24234422) 1 3x2 
ax (PrP zy P InP 


and analogous results for 0F)/dy and 0F3/0z. Adding the three derivatives results in 


3 3(x27+y?+27) 3 3 
.F = = _ 
VET BE FR re 


0, 


valid as long as |r| # 0, where F diverges. In electrostatics, F is proportional to the electric field of a 
point charge located at the origin. 
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Problems for Lecture 18 


1. Find the divergence of the following vector fields: 


a) F = xyi t yzj + zxk; 


b) F = yzi + xzj + xyk. 


Solutions to the Problems 


LECTURE 18. DIVERGENCE 


Lecture 19 | Curl 


View this lecture on YouTube 


Consider in Cartesian coordinates the three-dimensional vector field u = 41 (x,y,z)i + u2(x,y,z)j + 
uz(x,y,z)k. The curl of u, denoted as V x u and pronounced “del-cross-u”, is defined as the vector 


field given by 
B _ (Ms duy. (du dus\,., (duz ðu 
Vxu= 0/ox 0/dy 0/0z (s x ez a+ (SE xr. 
uy u2 u3 


Here, the cross product is used between a vector differential operator and a vector field. The curl 
measures how much a vector field rotates, or curls, around a point. A more math-based description 
will be given later. 

Example: Show that the curl of a gradient is zero, that is, V x (Vf) = 0. 

We have 


i j k 
Vx(Vf)=| 9/0x 09/9y 09/02 
af/ax af/ay af/dz 


LIT FT Vent OF cH ener. cf uon 
= (ayaz oz0y)  Vozox nz)?" VoxOy dyox) T| 


using the equality of mixed partials. 


Example: Show that the divergence of a curl is zero, that is, V - (V x u) — 0. 
We have 
_ ð (du dm\ 09 (du dus\ , A (duz du 
VAN x (5 SEHE dx) Oz \ ax oy 


m 9214 u un un l uz uz zi 
~ Noyoz 9z0y) | \dzdx  OxOz)]  NOxOy dydx) ” 


again using the equality of mixed partials. 
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Problems for Lecture 19 


1. Find the curl of the following vector fields: 


a) F = xyi t yzj + zxk; 


b) F = yzi + xzj dj 


- xyk. 


2. Consider a two-dimensional velocity field given by 


u = u(x, yji +u(x, y). 


Show that the vorticity w = V x u takes the form 


w=w(x, y)k. 


Determine w3 in terms of uz and up. 


Solutions to the Problems 


LECTURE 19. CURL 


Lecture 20 | Laplacian 


View this lecture on YouTube 


The Laplacian is the differential operator V - V = v2, given in Cartisian coordinates as 


The Laplacian can be applied to either a scalar field or a vector field. The Laplacian applied to a scalar 
field, f = f (x,y,z), can be written as the divergence of the gradient, that is, 


" | Of Pf Pf 
V?f — V.(Vf) alta wr 


The Laplacian applied to a vector field, acts on each component of the vector field separately. With 
u = u (x,y, zji + u»(x,y, z)3 + ua(x, y, z)k, we have 


V?u = V?uji + V?u»j + V?usk. 


The Laplacian appears in some classic partial differential equations. The Laplace equation, wave 
equation, and diffusion equation all contain the Laplacian and are given, respectively, by 


odo 


— 2 
3t DV*0. 


2 
vV? =0, = = gv, 


Example: Find the Laplacian of f (x,y,z) = x? +y* +2?. 


We have V?f 224-242 =6. 
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52 LECTURE 20. LAPLACIAN 
Problems for Lecture 20 


1 —— 
1. Compute V? C) forr #0. Here, r = yx? + y? +22. 


Solutions to the Problems 


Practice Quiz | The del operator 


1 
1. Let r = xi - yj - zk and r = yx? + y? +22. What is the value of V er 


2. Let F = =. The divergence V - F is equal to 


3 a IR 


An 


3. The curl of the position vector, V x r, is equal to 


r 

a) m 

b) 0 

T (V - r)r 
r 

d) i—-j+k 


Solutions to the Practice quiz 


53 


54 


PRACTICE QUIZ. THE DEL OPERATOR 


Lecture 21 | Vector derivative identities 


View this lecture on YouTube 


Let f = f(r) be a scalar field and u = u(r) and v = v(r) be vector fields, where u = ui + u»j + usk, 
and v = vi + vj + v3k. Here, we will change notation and define the position vector to be 
r = xyi c x53 + x3k (instead of using the coordinates x, y and z). We have already shown that 
the curl of a gradient is zero, and the divergence of a curl is zero, that is, 


V x Vf —0, V-(Vxu)=0. 
Other sometimes useful vector derivative identities include 
V. 


) 
) 
) 
V(u:v) 2 (u- V)v- (v: V)Jucux(Vxv)--vx(V xu), 
) 
) 


u:-V =u i Eu a u o 
1 Oxy | ?9x; | 33x3” 
which acts on a scalar field as 
m a Bi 
= idi n 9x1 "2 o ue Qx3' 


and acts on a vector field as 
(u- V)v = (u: Vr1)i - (u: Vv2) j + (wu: Voz) k. 


In some of these identities, the parentheses are optional when the expression has only one possible 
interpretation. For example, it is common to see (u: V)v written as u- Vv. The parentheses are 
mandatory when the expression can be interpreted in more than one way, for example V x u x v 
could mean either V x (u x v) or (V x u) x v, and these two expressions are usually not equal. 

Proof of all of these identities is most readily done by manipulating the Kronecker delta and Levi- 
Civita symbols, and I give an example in the next lecture. 
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LECTURE 21. VECTOR DERIVATIVE IDENTITIES 


Lecture 22 | Vector derivative identities 


(proof) 


View this lecture on YouTube 


To prove the vector derivative identities, we use component notation, the Einstein summation con- 
vention, the Levi-Civita symbol and the Kronecker delta. The ith component of the curl of a vector 


field is written using the Levi-Civita symbol as 


Quy 
(V xu); = EF 
and the divergence of a vector field is written as 


= Ou; 


We will continue to make use of the usual Kronecker delta and Levi-Civita symbol identities. 


As one example, I prove here the vector derivative identity 


V-(uxv)=v:-(Vxu)—u-(V xv). 


We have 
ð ; s . 
V.(uxv) = x (eikujor) (write using component notation) 

i 

Qu; dv 
= €ijk a Uy + EjjkUj 3x (product rule for the derivative) 

gie ) 
= veu; WEG Hei, Heer 
Dax TER, ij ip €ij ji 

=v-(Vxu)-u-(Vxv). (back to vector notation) 


The crucial step in the proof is the use of the product rule for the derivative. The rest of the proof just 


requires facility with the notation and the manipulation of the indices of the Levi-Civita symbol. 
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58 LECTURE 22. VECTOR DERIVATIVE IDENTITIES (PROOF) 


Problems for Lecture 22 


1. Use the Kronecker delta, the Levi-Civita symbol and the Einstein summation convention, and the 
identities 
a-b=0jajb;, (axb);—eiküjby, — ijk€ilm = Ójiókm — ÓjmÓki, 


to prove the following identities: 
a) V: (fu) -u- Vf - fV. u; 
b) Vx(Vxu)=V(V-u)- Vu. 
2. Consider the vector differential equation 


dr 
A = ult r(t), 


where 
T = xi X23 + x3k, u = ui + u53j + usk. 


a) Write down the differential equations for dx; /dt, dx2/dt and dx3/dt; 
b) Use the chain rule to determine formulas for d?x4 /dt?, d?x2/dt? and d*x3/dt?; 
c) Write your solution for d?r /dt? as a vector equation using the V differential operator. 


Solutions to the Problems 


Lecture 23 | Electromagnetic waves 


View this lecture on YouTube 


Maxwell's equations in free space are most simply written using the del operator, and are given by 


OE 
V.E=0, V.B=0, VxE=-, V x B = UE 


Here I use the SI units familiar to engineering students, where the constants ey and yo are called the 
permittivity and permeability of free space, respectively. 
From the four Maxwell's equations, we would like to obtain a single equation for the electric field 


E. To do so, we can make use of the curl of the curl identity 
V x(V x E) - VV. E) - VE. 


To obtain an equation for E, we take the curl of the third Maxwell's equation and commute the time 


and space derivatives 
9 
V x (V x E) = -4(V x B). 


We apply the curl of the curl identity to obtain 


9 
V(V-E)-V’E= —5,(V x B), 


and then apply the first Maxwell’s equation to the left-hand-side, and the fourth Maxwell’s equation 
to the right-hand-side. Rearranging terms, we obtain the three-dimensional wave equation given by 


QE 
38 = CVE, 


with c the wave speed given by c = 1/ ,/figéq ~ 3 x 10° m/s. This is, of course, the speed of light in 


Vacuum. 
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60 LECTURE 23. ELECTROMAGNETIC WAVES 


Problems for Lecture 23 


1. Derive the wave equation for the magnetic field B. 


Solutions to the Problems 


Practice Quiz | Vector calculus algebra 


1. Using the vector derivative identity V(u: v) = (u: V)v+(u-V)utux (V xv)+v x(V xu), 
determine which of the following identities is valid? 


a) P V(u:u) =u x (V x u) + (u: V)u 


b) ¿Y (uu) = u x (V x u) — (u: V)u 


c) 5 V(u-u) = (V xu) x u+ (u- V)u 


V(u:u)—(Vxu)xu-(u:V)u 
2. Which of the following expressions is not always zero? 
a) V x (Vf) 
b) V-(Vxu) 
o V-(Vf) 
d) V x (V(V -u)) 


3. Suppose the electric field is given by E(r,t) — sin(z — ct)i. Then which of the following is a valid 
free-space solution for the magnetic field B — B(r,t)? 


a) B(r,t) = * sin z — ct)i 
b) B(r,t) — ¿Sin z=ct)j 
c) B(r,t) = „sin x — ct)i 
d) B(r,t) — „sin x —ct)j 


Solutions to the Practice quiz 
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PRACTICE QUIZ. VECTOR CALCULUS ALGEBRA 


Week III 


Integration and Curvilinear 


Coordinates 


In this week's lectures, we learn about multidimensional integration. The important technique of using curvilin- 
ear coordinates, namely polar coordinates in two dimensions, and cylindrical and spherical coordinates in three 
dimensions, is used to simplify problems with circular, cylindrical or spherical symmetry. Differential operators 


in curvilinear coordinates are derived. The change of variables formula for multidimensional integrals using the 
Jacobian of the transformation is explained. 
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Lecture 24 | Double and triple integrals 


View this lecture on YouTube 


Double and triple integrals, written as 


[foam [f rosam | rav = fff, Fyz axayaz 


are the limits of the sums of AxAy (or AxAyAz) multiplied by the integrand. A single integral is the 


area under a curve y — f(x); a double integral is the volume under a surface z — f(x,y). A triple 
integral is used, for example, to find the mass of an object by integrating over its density. 

To perform a double or triple integral, the correct limits of the integral needs to be determined, 
and the integral is performed as two (or three) single integrals. For example, an integration over a 


rectangle in the x-y plane can be written as either 


yı x1 x yı 
/ ji f(xy)dxdy or / / f(x,y) dy dx. 
Vo Xo xo Vo 


In the first double integral, the x integration is done first (holding y fixed), and the y integral is done 
second. In the second double integral, the order of integration is reversed. Either order of integration 
will give the same result. 

Example: Compute the volume of the surface z — x?y above the x-y plane with base given by a unit square with 
vertices (0,0), (1,0), (1,1), and (0, 1). 

To find the volume, we integrate z — x?y over its base. The integral over the unit square is given by 


either of the double integrals 


1 pl 1 pl 
dE n x?y dx dy Or n / x?y dy dx. 
0 JO 0 JO 


The respective calculations are 


1 1 1 3 x=1 1 2 y=1 
n f Pydxdy = | ry dy = = ydy = 24 = 
o Jo o \ 3 lso 3 2 |, 
1 yl 1 2,/2 ¡y=1 1 fi 3¡x=1 
[ [omn (ae 
o Jo 0 y=0 0 a 


In this case, an even simpler integration method separates the x and y dependence and writes 


uL E 1 1 1 1 
eyaxdy= | ax | yay= (3) (5) =F 
ll M 0 quu 3J \2 6 
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66 LECTURE 24. DOUBLE AND TRIPLE INTEGRALS 


Problems for Lecture 24 


1. A cube has edge length L, with mass density increasing linearly from p1 to p» from one face of the 
cube to the opposite face. By solving a triple integral, compute the mass of the cube in terms of L, pı 


and p». 


Solutions to the Problems 


Lecture 25 Example: Double integral 


with triangle base 


View this lecture on YouTube 


Example: Compute the volume of the surface z = x*y above the x-y plane with base given by a right trian- 
gle with vertices (0,0), (1,0), (0,1). 


1 1 
0.8 0.8 
0.6 0.6 
> el 0 > y=1-rzr 
0.4 0.4 
0.2 0.2 
0 | 0 
0 0.2 04 06 08 1 0 0.2 04 06 08 1 
T "m 


The left figure illustrates the integral over x first and y second; and the right figure illustrates the 
integral over y first and x second. These are the two double integrals, 


1 pl-y 1 p1-x 
n / xy dx dy, or / | x?y dy dx. 
0 J0 0 JO 


The respective calculations are 


EY "(xy n 3 17 2 3,4 
L zydsdy- | E day=3 | Q-yydy=3 f W 3y +3y” — y") dy 
o Jo 0 x-ü 0 0 

1 


1 pa y 
312 "A 5 


1 pl-y 1 xy? y=1=x 1 A 1 fi 
n n Pydydx = | == da=5 | wand, | (x? — 2x + xt) dx 
0 v0 o \ 2 ly=0 2 Jo 2 Jo 
P. d 3v 
3 2 5) 


1 x3 E 1 
B M A. 60° 


x=1-y 


^ o1 
5-2 
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68 LECTURE 25. EXAMPLE: DOUBLE INTEGRAL WITH TRIANGLE BASE 


Problems for Lecture 25 


1. Compute the volume of the surface z — x?y above the x-y plane with base given by a parallelogram 
with vertices (0,0), (1,0), (4/3,1) and (1/3,1). 


Solutions to the Problems 


Practice Quiz | Multidimensional 


integration 


1. The volume of the surface z = xy above the x-y plane with base given by a unit square with vertices 
(0,0), (1,0), (1,1), and (0,1) is equal to 


NIe le Ale OF 


2. A cube has edge length of 1 cm, with mass density increasing linearly from 1g/ cm? to 2 g/ cm? 


from one face of the cube to the opposite face. The mass of the cube is given by 
a) 3.08 
b 15g 
c) 133g 
d) 1.0g 


3. The volume of the surface z — xy above the x-y plane with base given by the triangle with vertices 
(0,0), (1,1), and (2,0) is equal to 


wle Ale ORF om 


Solutions to the Practice quiz 
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PRACTICE QUIZ. MULTIDIMENSIONAL INTEGRATION 


Lecture 26 | Polar coordinates (gradient) 


View this lecture on YouTube 


In two dimensions, polar coordinates is the most commonly 
used curvilinear coordinate system. The relationship between 


Cartesian and polar coordinates is given by ó 
x —rcos6, y = rsin0. 
The unit vectors f and Ó are defined to be orthogonal and in 


the direction of increasing r and 0, respectively, and the position 


vector is given by r = rf. 


The f-Ó unit vectors are rotated an angle 0 from the i-j unit 


vectors. Trigonometry shows that 
f = cos 0i + sin0j, ô = —sin0i + cos 63. 


Notice that the direction of the unit vectors in polar coordinates are not fixed, but depend on their 


position. Here, £ = r(0) and 6 = Ó(0), and their derivatives are given by 


Now, consider a two-dimensional scalar field given by f — f(r,0). By using the definition of the 


two-dimensional gradient or by applying the chain rule, the differential df can be written as 


df=Vf-dr or df= Sarg Lao 


We have df 
dr = d(r?) = fdr + rao = fdr + rÓd0, 


so that equating the two forms for df results in 


V f - (fdr +r6d0) = Lar + F qe. 


By inspection, we see that 


or r 00 
so that the gradient is given by 
„ð ¿10 
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72 LECTURE 26. POLAR COORDINATES (GRADIENT) 


Problems for Lecture 26 


1. The inverse of a two-by-two matrix is given by 


Given 
f = cos i +sin0j, 6 = — sin 0i + cos 05, 


invert a two-by-two matrix to solve for ¿and j. 
2. Transform the Cartesian partial derivatives into polar form. 


a) Let f = f(x(r,0), y(r,0)). Using the chain rule and the definition of polar coordinates, show that 


of — 9f | . ,of of — wu m of 
3r coso | me 30 rsin-- | in 
b) Invert the result of Part (a) to show that 
u o3 sine ð T eunt , cos0 ð 
Ox arr 90 Oy rr 


3. Determine expressions for rf and r0 in Cartesian coordinates. 


Solutions to the Problems 


Lecture 27 | Polar coordinates 
(divergence & curl) 


View this lecture on YouTube 


Define a vector field in polar coordinates to be 
u = u;(r,0)f + ug(r,0)Ó. 


To compute the divergence and curl in polar coordinates, we will use # = f(0) and 6 = (0), and 


where k is the standard Cartesian unit vector pointing along the z-axis. 


The divergence is computed from 


IN ee 
Zn (^s. | 0795) u 


NE CMM AI A A 
=f mF (ug8) +0 735 f) $T. (ug8) 
| Ou, 0,1, ,109ug 

or Ex x 

129 , 190ug 

r or (ri) 4 r o0" 


" 210 
Vxu= c 5) x (up? + ug8) 
. 09 "P. ^ 109 FO 110 z 
=Px > f) +# x = (u00) +0 x (uf) +0 » 99 (00) 
Mar. A 10u.,. , 1 A 
— 04 AT (0x8) g (8xf)- us (8x) 


Note that the curl of a two-dimensional vector field defined in a plane points in the direction perpen- 
dicular to the plane following the right-hand rule. 
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74 LECTURE 27. POLAR COORDINATES (DIVERGENCE & CURL) 


Problems for Lecture 27 
1. Let u be a two-dimensional vector field given in polar coordinates by 


1 
u=- 
r 


(kır + k20) ; 


where kı and kə are constants. For r #0, determine V - u and V x u. 


Solutions to the Problems 


Lecture 28 | Polar coordinates 
(Laplacian) 


View this lecture on YouTube 


The two-dimensional Laplacian operator in polar coordinates can be computed using 
V*elvew), where V=r 


We have 


90 ¿109 90 ¿109 
2 a ‘ 

Vr = + @ HO 
(es 2) (es 2) 


aD 2 28 4 3X 313/3) 13/419 
"dr (+3) dr (6735) p^ y up (+3) er (6735) 


B^ uo: d gr. 1# 
9:22 ' rðr 799 


.183(,0|,19 
rorN ar) 12902 


The Laplacian operator can act on both a scalar or vector field. For the case of a vector field, one needs 


e 


to take into account the derivatives of the unit vectors with respect to 6. 
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76 LECTURE 28. POLAR COORDINATES (LAPLACIAN) 


Problems for Lecture 28 


1. Consider the fluid flow through a pipe of circular cross-section radius R with a constant pressure 
gradient along the pipe length. Define the z-axis to be the symmetry axis down the center of the pipe 
in the direction of the flowing fluid. The velocity field in a steady flow then takes the form u = u(r)k, 
where r = y/x2 + y? and u satisfies the Navier-Stokes equation given by 


Vu = a 
ur 

Here, G is the pressure gradient, v is the kinematic viscosity, and p is the fluid density, all assumed to 

be constant. You may further assume that the interior surface of the pipe has no slip so that the fluid 

velocity is zero when r = R. Solve for the velocity field u(r) in the pipe's cross section using the polar 

coordinate form for the Laplacian. What is the maximum value of the velocity? 


Solutions to the Problems 


Lecture 29 | Example: Central force 


View this lecture on YouTube 


A central force is a force acting on a point mass and pointing directly towards a fixed point in space. 
The origin of the coordinate system is chosen at this fixed point, and the axis orientated such that the 
initial position and velocity of the mass lies in the x-y plane. The subsequent motion of the mass is 
then two dimensional, and polar coordinates can be employed. 

The position vector of the point mass in polar coordinates is given by 


r — rf. 


The velocity of the point mass is obtained by differentiating r with respect to time. The added dif- 
ficulty here is that the unit vectors + = f(0(t)) and Ô = Ó(0(t)) are also functions of time. When 
differentiating, we will need to use the chain rule in the form 


dé dedo do, dÊ dódo de. 
de dodt dt’ dt dd d 


As is customary, here we will use the dot notation for the time derivative. For example, x = dx/dt 
and X = d?x/dt?. 
The velocity of the point mass is then given by 


CUP GET 
r= die = if 4- 7006; 
and the acceleration is given by 
dro „dô 
PEEL bree +100 + rd, 
= (f — r&?)e + (276 + r8)Ó 
A central force can be written as F = — ff, where usually f = f(r). Newton’s equation mi = F then 


becomes the two equations 
m(? — r6?) = — f, m(2i0 4- rà) = 0. 


The second equation is usually expressed as conservation of angular momentum, and after multipli- 


cation by r, is written in the form 


—(mr?ü) — 0, or mr? = constant. 
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78 LECTURE 29. CENTRAL FORCE 


Problems for Lecture 29 


1. The angular momentum I of a point mass m relative to an origin is defined as 
l=rxp, 
where r is the position vector of the mass and p = mr is the momentum of the mass. Show that 


[| = mr?|6]. 


dv d La ln 
vg a (2 


This result will be used to derive the conservation of energy. 


2. Prove that 


Solutions to the Problems 


Lecture 30 | Change of variables 


(single integral) 


View this lecture on YouTube 


A double (or triple) integral written in Cartesian coordinates, can sometimes be more easily computed 
by changing coordinate systems. To do so, we need to derive a change of variables formula. 
It is illuminating to first revisit the change-of-variables formula for single integrals. Consider the 


integral 
% 
i- f f (x) dx. 
Xo 


Let u(s) be a differentiable and invertible function. We can change variables in this integral by letting 
x = u(s) so that dx = u'(s) ds. The integral in the new variable s then becomes 


The key piece for us here is the transformation of the infinitesimal length dx = u’(s) ds. 
We can be more concrete by examining a specific transformation. Consider the calculation of the 
area of a circle of radius R, given by the integral 


R 
A-a[ V R2 — x? dx. 
0 


To more easily perform this integral, we can let x = R cos6 so that dx = —R sin 0 d0. The integral then 


becomes 


n/2 n/2 
A= ar? | y 1 — cos? O sin 0 dé = / sin? 0 d0, 
0 0 


which can be done using the double-angle formula to yield A — nR?. 
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LECTURE 30. CHANGE OF VARIABLES (SINGLE INTEGRAL) 


Lecture 31 | Change of variables 
(double integral) 


View this lecture on YouTube 


We consider the double integral 


I= COLS 


We would like to change variables from (x,y) to (s,t). For simplicity, we will write this change of 
variables as x = x(s,t) and y = y(s,t). The region A in the x-y domain transforms into a region A’ in 
the s-t domain, and the integrand becomes a function of the new variables s and t by the substitution 
f(x,y) = f(x(s,t),y(s,t)). We now consider how to transform the infinitesimal area dx dy. 

The transformation of dx dy is obtained by considering how an infinitesimal rectangle is trans- 
formed into an infinitesimal parallelogram, and how the area of the two are related by the absolute 
value of a determinant. The main result, which we do not derive here, is given by 


dx dy = |det (J)| ds dt, 


where J is called the Jacobian of the transformation, and is defined as 


d(x,y) _ [9x/ds 0x/0t 
dy/ds dy/0t)' 


To be more concrete, we again calculate the area of a circle. Here, using a two-dimensional integral, 


a= [f dx dy, 
A 


where the integral subscript A denotes the region in the x-y plane that defines the circle. To perform 


the area of a circle can be written as 


this integral, we can change from Cartesian to polar coordinates. Let 


x=rc0s0, y = rsin0. 
We have 
ax/dr 0x/90 0 —rsin0 
ee is al Aa codd 
dy/dr 0y/90 sin 0 T cos 0 


The region in the r- plane that defines the circle is 0 € r < Rand 0 € 0 < 27. The integral then 


becomes 
270 R 27 R 
A= | n rardo = | ae | rdr = mR. 
0 0 0 0 
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82 LECTURE 31. CHANGE OF VARIABLES (DOUBLE INTEGRAL) 


Problems for Lecture 31 


1. The mass density of a flat object can be specified by 7 = c(x, y), with units of mass per unit area. 
The total mass of the object is found from the double integral 


M= | | ot.waray. 


Suppose a circular disk of radius R has mass density po at its center and p at its edge, and its density 
is a linear function of the distance from the center. Find the total mass of the disk. 


00 


2. Compute the Gaussian integral given by I = / e". dx. Use the well-known trick 


— 00 


00 2 co co co 00 
P= (/ e dx) =f e? dx | e Y ay = | f e +9") dx dy. 


Solutions to the Problems 


Practice Quiz | Polar coordinates 


1. rÓ is equal to 
a) xi-- yj 
b) xi — yj 
c) yi t xj 
d) —yi + xj 


A 


2 a lt 
- jg Y equal to 


3. Suppose a circular disk of radius 1 cm has mass density 10 g/cm? at its center, and 1 g/cm? at its 
edge, and its density is a linear function of the distance from the center. The total mass of the disk is 


equal to 
a) 8.80 g 
b) 10.21 g 
c) 12.57 g 
d) 1723 g 


Solutions to the Practice quiz 
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PRACTICE QUIZ. POLAR COORDINATES 


Lecture 32 | Cylindrical coordinates 


View this lecture on YouTube 


Cylindrical coordinates extends polar coordinates to three 
dimensions by adding a Cartesian coordinate along the z- 
axis. To conform to standard usage, we change notation 
and define the relationship between the Cartesian and the 
cylindrical coordinates to be 


x —pcosq, y = psin $, Z=2Z; 
with inverse relation 
p= V3 y), tang = y/x. 


A spatial point (x,y,z) in Cartesian coordinates is now 


specified by (p, Q, z) in cylindrical coordinates (see figure). 


The orthogonal unit vectors o, $, and 2 point in the direction of increasing p, $ and z, respectively. 
The position vector is given by r = pp + z2. The volume element transforms as dx dy dz = p dp de dz. 

The del operator can be found using the polar form of the Cartesian derivatives (see the problems). 
The result is 


The Laplacian, V - V, is computed taking care to differentiate the unit vectors with respect to ¢: 


v? 183,195 0 
9p pop pag? 922 


4337.9. 1# @ 
pap V 3p) "ap? ag 


The divergence and curl of a vector field, A = Ap(p, P,2)P + Ag(p, $,z)ó + Az(p, $, z)2, are given by 
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86 LECTURE 32. CYLINDRICAL COORDINATES 


Problems for Lecture 32 


1. Determine the del operator V in cylindrical coordinates. There are several ways to do this, but a 


straightforward, though algebraically lengthy one, is to transform from Cartesian coordinates using 


$ = cos pô — sin pd, Y = sin gp + cos pd, 


and 
sing ð 


p 9$ 


cos $ ð 
p op 


pm T 
y dy N 


d ð 
a cos $ ap 

2. Compute V - 6 in two ways: 
a) With the divergence in cylindrical coordinates; 

b) By transforming to Cartesian coordinates. 


3. Using cylindrical coordinates, compute V x 6, V - $ and V x @. 


4. The center-of-mass of a solid with density p and total mass M is defined (with respect to a given 


coordinate system) as 


1 
== dV. 
úl s f," 


Find the center-of-mass of the uniform solid cone pictured below, with coordinate system specified. 


1 
You may assume that the volume of the cone is given by V — gra. 


Solutions to the Problems 


Lecture 33 | Spherical coordinates 
(Part A) 


View this lecture on YouTube 


Spherical coordinates are useful for problems with 
spherical symmetry. The relationship between the Carte- 


sian and the spherical coordinates is given by 


x —rsin6cos ¢, y = rsinOsinf, z —rcosÓ, 


where 0 < 0 < m and 0 € $ < 27m. A spatial point 
(x,y,z) in Cartesian coordinates is now specified by 
(r,0, $) in spherical coordinates (see figure). The orthog- 
onal unit vectors f, Ô, and $ point in the direction of 


increasing r, 0 and @, respectively. 

The position vector is given by 

and the volume element transforms as 
dx dy dz = r° sin 0 drd0 dọ. 

The spherical coordinate unit vectors can be written in terms of the Cartesian unit vectors by 


f = sin 0 cos $ à + sin 0 sin $ j + cos 0 k, 
= cos cos $ à + cos O sinp j — sin 0 k, 


ô 
$ — —sin $i + cos $ j; 
with inverse relation 


i = sin cos p f + cos0 cos p Ô — sing à, 
j = sin O sin $ f + cos O sin $ Ó + cos pd, 


k = cos 0 f — sind Ô. 


By differentiating the unit vectors, we can derive the sometimes useful identities 


oF, 96 A d$ 
eo gh a” 
P A 96 " 909 — A 
T = sin 0$, d = cos 06, 36 = — sin 0f — cos 00. 
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88 LECTURE 33. SPHERICAL COORDINATES (PART A) 


Problems for Lecture 33 


1. Write the relationship between the spherical coordinate unit vectors and the Cartesian unit vectors 
in matrix form. Notice that the transformation matrix Q is orthogonal and invert the relationship 
using Q7! = QF. 


2. Use the Jacobian change-of-variables formula for triple integrals, given by 


ox/dr 0x/900 dx/dp 
dx dy dz = |det | dy/or 9y/96 0y/0p || dr dð dp, 
Oz/ór 02/08 0z/0$ 


to derive dx dy dz = 1? sin 0 dr d0 dọ. 


3. Consider a scalar field f = f(r) that depends only on the distance from the origin. Using dx dy dz = 
1? sin 0 dr d0 dp, and an integration region V inside a sphere of radius R centered at the origin, show 
that 
R 
/ fdV = an | r f(r) dr. 
V 0 

4. Suppose a sphere of radius R has mass density pọ at its center, and p at its surface, and its density 
is a linear function of the distance from the center. Find the total mass of the sphere. What is the 
average density of the sphere? 


Solutions to the Problems 


Lecture 34 | Spherical coordinates 
(Part B) 


View this lecture on YouTube 


First, we determine the gradient in spherical coordinates. Consider the scalar field f = f(r,0,4). 
Our definition of a total differential is 


_ of, of Usa, 
df = ydri 5540 + = Vf dr. 


In spherical coordinates, 


and using 
FL or P 
—=0 =~ = sin0 
Ti ag u) > 
we have 
ESEL dd. = drê 4-rd0Ó - rsinüdo do 
CS BIET TRES Td f sind de ©. 


Using the orthonormality of the unit vectors, we can therefore write df as 


f. Lat, 1 21, EU A R 
df = E 79997 | pane oe! (dr? -- rd8 0 -- rsin6d$ 9), 


showing that 
vi-29f5,19fg, 1 of 


or | roo reine og 


Some messy algebra will yield for the Laplacian 


10 / ðf 1 9f 1 æf 
2 

B Ä eZ | 4 
Via y (A L) T2 sin 90 s (sin 5) r2 sin? 0 99? 


and for the divergence and curl of a vector field, A = Ay(r,0,p)? + Ae(r,0,$)0 + Ag(r,0,$)®, 


BICI 1 04, 
dS gar r) maa PP A Tolg y dp ' 
f 9 9A] Ó[ 1 OA, $ 9A, 
A= Ay) rA) | 4 A 
A nS ES 9) 1 r E ap x o) r E (r46) | 
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90 LECTURE 34. SPHERICAL COORDINATES (PART B) 


Problems for Lecture 34 


1. Using the formulas for the spherical coordinate unit vectors in terms of the Cartesian unit vectors, 


prove that 
or x of . A 
38 = 0, Em = sin Q. 


2. Compute the divergence and curl of the spherical coordinate unit vectors. 


3. Using spherical coordinates, show that for r + 0, 


eu 


Solutions to the Problems 


Practice Quiz Cylindrical and 
spherical coordinates 


1. With p = J/x? 4-32, the function V? (5) is equal to 


2. Let r = xi. Then (%, Ó, à) is equal to 
a) (i,j,k) 
b) (i, k, 5) 


c) i,—k,j) 


d) i, k, —j) 


3. Suppose a sphere of radius 5cm has mass density 10 g/cm? at its center, and 5g /cm? at its surface, 
and its density is a linear function of the distance from the center. The total mass of the sphere is given 


by 
a) 3927 g 
b) 3491 g 
c) 3272g 
d) 3142g 


Solutions to the Practice quiz 
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PRACTICE QUIZ. CYLINDRICAL AND SPHERICAL COORDINATES 


Week IV 


Line and surface integrals 


This week we learn about line integrals and surface integrals. We learn how to take the line integral of a scalar 
field and use line integrals to compute arc lengths. We then learn how to take line integrals of vector fields by 
taking the dot product of the vector field with tangent unit vectors to the curve. Consideration of the line integral 
of a force field results in the work-energy theorem. Next, we learn how to take the surface integral of a scalar field, 
and compute the surface areas of a cylinder, cone, sphere and paraboloid. We then learn how to take the surface 
integral of a vector field by taking the dot product of the vector field with normal unit vectors to the surface. The 
surface integral of a velocity field is used to define the mass flux of a fluid through the surface. 
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Lecture 35 | Line integral of a scalar field 


View this lecture on YouTube 


Define the line integral of a scalar field f = f(r) over a curve C by subdividing the curve into 
small scalar elements of length ds, multiplying each ds by the average value of f on the element, and 


summing over all elements. We write the line integral as 


OL 


For simplicity, we restrict our discussion here to curves in the x-y plane. By the Pythagorean theorem, 
we can write 

ds = \/ (dx)? + (dy), 
and it is possible to convert a line integral into an ordinary one-dimensional integral in one of two 
ways. First, if the curve can be specified by a one-dimensional function, such as y = y(x) from x = xo 
to xr, then dy = y'(x)dx and 

ds = 1/1+y'(x) dx. 


] fes - [^ f(x, y(x ) V1 ey» y dx. 


Second, if the curve can be parameterized by t, such as x = x(t) and y = y(t), with t ranging from 
t = to to tr, then dx = x(t)dt and dy = y(t)dt, and 


We then have 


ds = x(t? + y(t)? dt. 


The line integral then becomes 


[ f(x,y) ds = MEO x)? + y (0? dt. 


95 


96 LECTURE 35. LINE INTEGRAL OF A SCALAR FIELD 


Problems for Lecture 35 


1. Suppose that a curve can be specified by the one-dimensional function y — y(x). Parameterize the 
curve to show that the two forms of the line integral are equivalent. 


Solutions to the Problems 


Lecture 36 | Arc length 


View this lecture on YouTube 


The arc length P of a curve C is given by 


P= f ds 
C 


For even simple curves, this integral may not have an analytical solution. For example, consider the 
circumference of an ellipse. The equation for an ellipse centered at the origin with width 2a and height 
2b is given by 


2 2 
Cal 


= 
Assuming a > b, the eccentricity of this ellipse is defined as 


x 
at 
pes yı — (b/a)2. 
To compute the circumference of an ellipse, we parameterize the ellipse using 
x(0) = acos6, y(0) = bsin0, 


where 6 goes from 0 to 27t. The infinitesimal arc length ds is then given by 


ds = \/ (dx)? + (dy)? = \/x’(0)2 + y (0)? do 
= Va? sin? 0 + P? cos? 0 d0 


= aA — [1 — (b/a)?] cos? 6 d8 
=av1 — e? cos? 0 d0. 


Therefore, the circumference of an ellipse — or perimeter P — is given by the line integral 


270 
p= fds=af v 1 — e cos? 040 
C 0 
m/2 


= 4a V1 — e? cos? 040, 
0 


where the final integral is over one-quarter of the arc length of the ellipse , and is called the complete 
elliptic integral of the second kind. 
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98 LECTURE 36. ARC LENGTH 


Problems for Lecture 36 


1. Compute the circumference of a circle of radius R by parameterizing the circle and performing a 
line integral. 


Solutions to the Problems 


Lecture 37 | Line integral of a vector field 


View this lecture on YouTube 


Let u(r) be a vector field, C a directed curve, and dr the infinitesimal displacement vector along 
C. Define a unit vector f such that dr = fds. Then the line integral of u along C is defined to be 


[cara | w tas, 
E C 


which is a line integral of the scalar field u - t. If the curve is closed, sometimes a circle is written in 
the middle of the integral sign. 

A general method to calculate the line integral of a vector field is to parameterize the curve. Let 
the curve be parameterized by the function r = r(t) as t goes from to to ty. Using dr = (dr /dt)dt, the 


[ua = Puro) i dt. 


to 


line integral becomes 


Sometimes the curve is simple enough that dr can be computed directly. 


Example: Compute the line integral of r = xà + yj in the x-y plane along two curves from the origin to the 
point (x,y) — (1,1). The first curve C4 consists of two line segments, the first from (0,0) to (1,0), and the 
second from (1,0) to (1,1). The second curve C» is a straight line from the origin to (1,1). 


The computation along the first curve requires two separate integrations. For the curve along the 


x-axis, we use dr = dxi and for the curve at x = 1 in the direction of j, we use dr = dyj. The line 


1 1 
[vara fm f ydy — 1. 
Ci 0 0 


For the second curve, we parameterize the line by r(t) = t(i +j) as t goes from 0 to 1, so that 


integral is therefore given by 


dr = dt(i + j), and the integral becomes 


1 
f rdr = f 2tdt = 1. 
C5 0 


The two line integrals are equal, and for this case only depend on the starting and ending points of 
the curve. 


99 


100 LECTURE 37. LINE INTEGRAL OF A VECTOR FIELD 


Problems for Lecture 37 


1. In the x-y plane, calculate the line integral of the vector field u — —yi + xj counterclockwise around 
a square with vertices (0,0), (L,0), (L, L), and (0, L). 


2. In the x-y plane, calculate the line integral of the vector field u = —yö + xj counterclockwise around 


the circle of radius R centered at the origin. 


Solutions to the Problems 


Lecture 38 | Work-energy theorem 


View this lecture on YouTube 


Newton's second law of motion for a mass m with velocity v acted on by a force F is given by 


The conservation of energy is an important concept in physics and requires a definition of work. We 
first take the dot product of both sides of Newton's law with the velocity vector and use the identity 


to obtain 


Integrating from an initial time t; to a final time ty, with v(t;) = v; and v(ty) = vf, we have 


1 ; 1 N tf 
lve mios] =| P-vdt= | Far, 


i 


where we have used dr — vdt. The line integral is taken along the curve C traversed by the mass 


between the times t; and t f. We define the kinetic energy T of a mass m by 
1 
T — 2m 2 
"lel 
and the work W done by a force on a mass as it moves along a curve C as 


w= | Fear, 
C 


We these definitions, the work-energy theorem states that the work done on a mass by a force is equal 
to the change in the kinetic energy of the mass, or in equation form, 


W =T; - T;. 
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102 LECTURE 38. WORK-ENERGY THEOREM 


Problems for Lecture 38 


1. A mass m is dropped from a building of height h. Assuming the mass falls with constant accel- 
eration —$, calculate the work done by gravity and compute the velocity of the mass as it hits the 


ground. 


Solutions to the Problems 


Practice Quiz | Line integrals 


1. The arc length of the parabola y = x? for 0 < x < 1 is given by the integral 


a) [ V1+2x dx 


b) [ V1+4x dx 


1 
c) / V1+2x2 dx 
0 
1 


d) I V1+4x2 dx 


2. The line integral of the vector field u = —yi + xj counterclockwise around a triangle with vertices 
(0,0), (L,0), and (0, L) is equal to 


d) 212 


3. A mass m is shot upward from the ground with a speed vo, attains a maximum height, and then 
falls back to ground. Calculate the work done by gravity on the mass. 


a) 0 

b) 05/48 
c) 05/28 
d) 05/8 


Solutions to the Practice quiz 
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104 PRACTICE QUIZ. LINE INTEGRALS 


Lecture 39 | Surface integral of a scalar 
field 


View this lecture on YouTube 


Define the surface integral of a scalar field f = f(r) over a surface S by subdividing the surface into 
small elements d5, multiplying each element by the average value of f on the element, and summing 
over all elements. We write the surface integral as 


[ f(r) dS. 


A parameterization of the surface S is given by 
r(u,v) = x(u,v)i +y(u,v)j 4 z(u,v)k, 


where the position vector r points to the surface, and the surface is spanned as the parameters u and 
v vary. The surface integral is performed over u and v, and the infinitesimal surface element dS is 
found from the area defined by the tangent vectors dr/du and dr/dv to the surface by way of the 


cross product, that is, 


or or 
— X — 


gue ðu dv 


du dv. 


If the surface lies over the x-y plane and can be describe by z = z(x, y), then we can use x and y as our 
parameters and write 
r(x,y) = xi 4- yj t z(x, y)k. 


Therefore, 


Seren u AR N 
dy 


and the cross product is given by 


e edo o az. ə 
= es 0 dz/ax| = -Ži Zi+k 
Ox Oy oy 

0 1 0z/0y 


The surface integral in this case becomes 


[res f [reis y) | (=) | (5) iav 
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106 LECTURE 39. SURFACE INTEGRAL OF A SCALAR FIELD 


Problems for Lecture 39 


1. Compute the lateral surface area A = f, dS of a cylinder (see figure) in two ways. 


4) Unroll the cylinder and compute the area of the resulting rectangle. 


b) Define the cylinder parametrically as 
r =acos0i+asin0j+zk, forO<z<b and 0<@<2z7, 


and compute the surface integral. 


2. Compute the lateral surface area A = f dS of a cone (see figure) in two ways. 


a) Unroll the cone and compute the area of the resulting circular sector. 


b) Define the cone parametrically as 


r= " costi - = sind j +zk, forO<z<b and 0<0<2T, 


and compute the surface integral. 


Solutions to the Problems 


Lecture 40 | Surface area of a sphere 


View this lecture on YouTube 


Surface area is given by the integral 


s- | ds, 
S 


where the integral is evaluated by parameterizing the surface. Here, we compute the surface area of a 
sphere of radius R. The sphere can be parameterized using spherical coordinates as 


r(0,p) = Rsin0cos pi + R sin 8 sin pj + RcosOk, for 0<0<T1 and O<d<2n. 


To find the infinitesimal surface element, we compute the partial derivatives of r: 


d 
T = Rcos8cos i+ Rcos0sin $ j — Rsin6 k, 9i = —Rsin f sin $$ + Rsin 0 cos ¢ j. 
The cross product is 
: j k 
or or ? : . . 
59 * 20 —|Rcos0cosp  Rcos0sin$ —Rsin6O 
—Rsin@sing Rsindcos$® 0 
= R?sin? 0 cos p i + R? sin? 0 sin p j + R2 sin cos k, 
so that 
ar x cu NE R? (sint 8 cos? p + sin! 8 sin? y + sin? 6 cos? 8 
90 0$ 
= R?sin0 V sin? 0(sin? $ + cos? p) + cos? 0 
= R?sin6. 
Therefore, 
or gr 2. 
dS — 36 * ap d0 d$ = R^ sin 0 d0 d. 


Note that for a sphere, dV = 1? sin @ dr d0 dọ, so that dSdr is equal to dV at r = R. 
The surface area of a sphere can be found from the integral 


2m pr 27 T 
s- [as- f n Rsinododp=R | dp f sin 0 d0 = Ax R?. 
S 0 0 0 0 
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108 LECTURE 40. SURFACE AREA OF A SPHERE 
Problems for Lecture 40 
1. Compute the lateral surface area S = [, dS of a paraboloid (see figure), defined by 
wz =b (x? + 2) for0<z<b 
= y), or0<z<b. 


Compute the resulting two-dimensional integral using polar coordinates. 


Solutions to the Problems 


Lecture 41 | Surface integral of a vector 
field 


View this lecture on YouTube 


Let u = u(r) be a vector field, and let S be a surface with normal unit vector f. The surface in- 


fuas = [u nas, 
S 5 


which is a surface integral of the scalar field u f». The integration surface can be either open or closed. 


tegral of u over S is defined to be 


For an open surface, the direction of the normal vectors needs to be specified (such as up or down), 
but for a closed surface, n is always assumed to be in the outward direction. 
Example: Compute the surface integral of r = xà + yj + zk over a cube centered at the origin with sides parallel 


to the axes and side lengths equal to L . 


The faces of the cube are located at x = +L/2, y = +L/2 and z = +L/2, and the surface area of each 


face is L?. For the face at x = L/2, say, we have r = (L/2)i + yj +zk and dS = idS. The integral over 


L L? 
rds => fds ==. 
[ 2 Js 2 


One can see that the surface integrals over all six faces of the cube are equal, and we obtain for the 


this face is given by 


surface integral over the entire cube, 
L3 
pr-ds =6x Zul, 
S 2 


equal to three times the volume of the cube. The circle notation on the integral sign signifies integration 


over a closed surface. 


Example: Compute the surface integral of r = xà + yj + zk over a sphere centered at the origin with radius R . 


Using spherical coordinates, on the surface of a sphere of radius R centered at the origin, we have 
r = Rf, dS = fdS and the surface area of the sphere is Az R?. Therefore, 


rias Ras = ante, 
S S 


equal to three times the volume of the sphere. 
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110 LECTURE 41. SURFACE INTEGRAL OF A VECTOR FIELD 


Problems for Lecture 41 


1. Compute the surface integral of r = xi + yj + zk over a closed cylinder centered at the origin with 


with radius a and length 1. 


Solutions to the Problems 


Lecture 42 | Flux integrals 


View this lecture on YouTube 


The surface integral of a vector field is often called a flux integral. If u is the fluid velocity (length 
divided by time), and p is the fluid density (mass divided by volume), then the surface integral 


[owas 
S 


computes the mass flux, that is, the mass passing through the surface 5 per unit time. Flux integrals 


for the electric and magnetic vector fields are also defined. 
Example: Find the flux of the electric field through a sphere of radius R centered at the origin, where a point 
charge q is located. The electric field due to the point charge is given in spherical coordinates by Coulomb's law, 


q 


- f. 
Artegr2 


To compute the flux integral, we use dS = fdS, where the surface area of the sphere is given by 


S = 471R2. We have 
q q 
E.dS = dS = 
$ = m? ? €0” 


which is observed to be independent of R because of the Coulomb inverse square law. 
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112 LECTURE 42. FLUX INTEGRALS 


Problems for Lecture 42 


1. Calculate the mass flux of a laminar fluid of density p, viscosity v and constant pressure gradient 
G passing through a cross section of a pipe of radius R. Choosing z as the symmetry axis for the pipe, 
the velocity of the fluid is given by 


u(r) =m (1- (2) a 


where r is the radial coordinate in the cross section and 


_ GR? 
~ 4vp 


Um 


is the maximum velocity of the fluid in the center of the pipe. 


Solutions to the Problems 


Practice Quiz | Surface integrals 


1. The shape of a donut or bagel is called a torus. Let R be the radius from the center of the hole to the 
center of the torus tube and r be the radius of the torus tube. Then the equation for a torus symmetric 


e 
(r yer) rz =p, 


The torus may be parameterized by 


about the z-axis is given by 


x = (R4 rcos0) cos ọ, y = (R+rc0s0) sing, z —rsin6, 
where 0 < 0 < 27 and 0 < $ < 27. The infinitesimal surface element dS for the torus is given by 
a) R(r + R cos q) d0 dp 
b) r(R+rcos q) dé d 
c) R(r + Rcos0) d0 d 
d) r(R +r cos 8) de d 


2. Consider a right circular cylinder of radius R and length L centered on the z-axis. The surface 


integral of u = xi + yj over the cylinder is given by 
a) 0 
b) TR?L 
c) 2x R?L 
d) Ax R?L 


3. The flux integral of u = zk over the upper hemisphere of a sphere of radius R centered at the origin 


with normal vector f is given by 


d) AR? 


Solutions to the Practice quiz 
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114 PRACTICE QUIZ. SURFACE INTEGRALS 


Week V 


Fundamental Theorems 


In this week's lectures, we learn about the fundamental theorems of vector calculus. These include the gradient 
theorem, the divergence theorem, and Stokes' theorem. We show how these theorems are used to derive the law of 


conservation of energy, continuity equations, define the divergence and curl in coordinate-free form, and convert 
the integral version of Maxwell's equations to differential form. 
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Lecture 43 | Gradient theorem 


View this lecture on YouTube 


The gradient theorem is a generalization of the fundamental theorem of calculus for line integrals. 
Let V $ be the gradient of a scalar field p = $(r), and let C be a directed curve that begins at the 
point rı and ends at r2. Suppose we can parameterize the curve C by r = r(t), where tı < t < to, 
r(t) = ri, and r(t?) = r2. Then using the chain rule in the form 


d dr 
Solr) = vete) Sr, 


and the standard fundamental theorem of calculus, we have 


b dr 2 q 
| ve - f ve) maro ee 


= $(r(t3)) — p(r(t1)) = Pr) - (n). 
A more direct way to derive this result is to write the differential 
dp = Ve - dr, 


so that 
[vo-ar = | ap = ola) - or). 
C E 


We have thus shown that the line integral of the gradient of a function is path independent, depending 
only on the endpoints of the curve. In particular, we have the general result that 


f vean co 


for any closed curve C. 
Example: Compute the line integral of r = xà + yj in the x-y plane from the origin to the point (1,1). 
We have r = 1 V (x? + y?), so that the line integral is path independent. Therefore, 

(11) 


[rar =; vem ee =1. 
c 2 Jc 2 (0,0) 
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118 LECTURE 43. GRADIENT THEOREM 


Problems for Lecture 43 
1. Let (r) = x?y + xy? +2. 
a) Compute VE. 
b) Compute fe Vp dr from (0,0,0) to (1,1, 1) using the gradient theorem. 
c) Compute f- Vo -dr along the lines segments (0,0,0) to (1,0,0) to (1,1,0) to (1,1, 1). 


Solutions to the Problems 


Lecture 44 | Conservative vector fields 


View this lecture on YouTube 


For a vector field u defined on RÌ, except perhaps at isolated singularities, the following conditions 


are equivalent: 
1. Vxu=0; 
2. u = V $ for some scalar field $ = p(r); 


3. f u-dr is path independent for any curve C; 
C 


4. $ u- dr = 0 for any closed curve C. 
C 


When these conditions hold, we say that u is a conservative vector field. 


Example: Let u(x,y) = x?(1 + y?)i + y? (1 + x?)3. Show that u is a conservative vector field, and determine 
$ = $(x,y) such that u = V4. 
To show that u is a conservative vector field, we can prove V x u = 0: 
i j k 
Vxu=| 9/9x 9/dy 9/0 = (ey? — 3x74? )k =0. 
lv) FO) 0 


To find the scalar field $, we solve 


a =x7(1+y°), x =y (1417). 


Integrating the first equation with respect to x holding y fixed, we find 


p= Parra $0 9) fo), 


where f = f(y) is a function that depends only on y. Differentiating $ with respect to y and using the 
second equation, we obtain 


xyv.f())-yü0-xx3) o  f(y-y. 


One more integration results in f (y) = 1?/3 +c, with c constant, and the scalar field is given by 


1 
$6, y) = (Ptr ty’) e. 
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120 


Problems for Lecture 44 


1. Let u = (2xy +22) + (yz + x?)3 + (2zx 4 


a) Show that u is a conservative vector field. 


b) Calculate the scalar field $ such that u = V $. 


Solutions to the Problems 


LECTURE 44. CONSERVATIVE VECTOR FIELDS 


Lecture 45 | Conservation of energy 


View this lecture on YouTube 


The work-energy theorem states that the work done on a mass by a force is equal to the change 
in the kinetic energy of the mass, or 


Lre-5- 
where the kinetic energy of a mass m moving with velocity v is given by 
T = Emol. 
If F is a conservative vector field, then we can write 
F=-VV, 
where V = V(r) is called the potential energy. Using the gradient theorem, we have 
-h= | VV-dr=Vi— Vy, 
where V; and V; are the initial and final potential energies of the mass. Rearranging terms, we have 
T; + Vi = Ty + Vy. 


In other words, the sum of the kinetic and potential energy is conserved. 
Example: Find the potential energy of a mass m in the gravitational field of a mass M. 


We place the origin of our coordinate system on the mass M. The gravitational force on m at position 


r is then given by the inverse square law, written as 


mMr 


F=-G 
y3 


In the problems of Lecture 17, we have shown that V(1/r) = —r/ 1?. Therefore, the potential energy 


of m is given by 
V= ud 
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122 LECTURE 45. CONSERVATION OF ENERGY 


Problems for Lecture 45 


1. The escape velocity is the smallest initial velocity for a mass on the Earth's surface to escape from 
the Earth's gravitational field. Using the conservation of energy, determine the escape velocity of a 
mass m. Define M to be the mass of the Earth, and R its radius. The gravitational constant G and the 
acceleration due to gravity on the surface of the Earth g are related by 


_ GM 
8= RE 


Write the escape velocity in terms of g and R. 


Solutions to the Problems 


Practice Quiz | Gradient theorem 


1. Let $(r) = xyz. The value of | V $ - dr from (0,0,0) to (1,1, 1) is equal to 
C 


a) (x+y)? +z 


Fxy ty +z 


d) x? — xy +y? +z 


Solutions to the Practice quiz 
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124 PRACTICE QUIZ. GRADIENT THEOREM 


Lecture 46 Divergence theorem 


View this lecture on YouTube 


Let u be a differentiable vector field defined inside and on a smooth closed surface S enclosing a 
volume V. The divergence theorem states that the integral of the divergence of u over the enclosed 


volume is equal to the flux of u through the bounding surface; that is, 


| 70 dv =$ uds. 


We first prove the divergence theorem for a rectangular solid with sides parallel to the axes. Let 
the rectangular solid be defined bya < x < b,c < y € d, and e € z < f. With u = ui + u2j + usk, 


the volume integral over V becomes 


je u) dV = [ff a eme) dx dy dz. 


The three terms in the integral can be integrated separately using the fundamental theorem of calculus. 


Each term in succession is integrated as 


F [ ne P dy J dyaz= f fi (ur (b, y, z) — ux(a, y, z)) dy dz; 


F [ (f 7 a) iis - f f (u»(x,d,z) — u»(x,c,z)) dx dz; 
pu (f a) dxdy = [f (ua (x, y, f) — us(x,y,e)) dx dy. 


The integrals on the right-hand-sides correspond exactly to flux integrals over opposite sides of the 
rectangular solid. For example, the side located at x = b corresponds with dS = i dy dz and the side 
located at x = a corresponds with dS = —idy dz. Summing all three integrals yields the divergence 
theorem for the rectangular solid. 

Now, given any volume enclosed by a smooth surface, we can subdivide the volume by a very fine 
three-dimensional rectangular grid and apply the above result to each rectangular solid in the grid. All 
the volume integrals over the rectangular solids add. The internal rectangular solids, however, share 
connecting side faces through which the flux integrals cancel, and the only flux integrals that remain 
are those from the rectangular solids on the boundary of the volume with outward facing surfaces. 


The result is the divergence theorem for any volume V enclosed by a smooth surface 5. 
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126 LECTURE 46. DIVERGENCE THEOREM 


Problems for Lecture 46 


1. Prove the divergence theorem for a sphere of radius R centered at the origin. Use spherical coordi- 


nates. 


Solutions to the Problems 


Lecture 47 | Divergence theorem 


(example 1) 


View this lecture on YouTube 


The divergence theorem is given by 


| v iav = ducas. 


Test the divergence theorem using u = xy -F- yz j + zx k for a cube of side L lying in the first octant with a 


vertex at the origin. 


Here, Cartesian coordinates are appropriate and we use V - u = y +z + x. We have for the left-hand 


side of the divergence theorem, 


[vmm [p orem: 


=1*/241*/241*/2 
=31*/2. 
For the right-hand side of the divergence theorem, the flux integral only has nonzero contributions 


from the three sides located at x = L, y = Land z = L. The corresponding unit normal vectors are i, 


j and k, and the corresponding integrals are 


L pL L pL L pL 
puas= f J Lydydz+ | n Lzdxdz | f Lx dx dy 
S 0 JO 0 JO 0 J0 


= 15/2 + 13/2 + 13/2 
= 3154/2. 
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128 LECTURE 47. DIVERGENCE THEOREM (EXAMPLE 1) 


Problems for Lecture 47 


1. Test the divergence theorem using u = x?y ¿+ y?z j +22x k for a cube of side L lying in the first 


octant with a vertex at the origin. 


2. Compute the flux integral of r = xi + yj + zk over a square box with side lengths equal to L by 
applying the divergence theorem to convert the flux integral into a volume integral. 


Solutions to the Problems 


Lecture 48 Divergence theorem 


(example 2) 


View this lecture on YouTube 


The divergence theorem is given by 


] v iav = ducas. 


Test the divergence theorem using u = r?f for a sphere of radius R centered at the origin. 


To compute the left-hand-side of the divergence theorem, we recall the formula for the divergence of 


a vector field u in spherical coordinates: 


1 oO 1 dup 


10 : 
(sin Hug) + rend ap š 


P 2 
va art ur) 4 rsin 0 90 


Here, u, — r? is the only nonzero component of u, and we have 


Vu= 54 (r) — 4r. 


Therefore, using dV — r? sin 0 dr d0 dq, we have 


27 pm pR 
[| va f / / 4r? sin 0 dr d0 dp 
V 0 0 J0 


27 T R 
= ae f sinodo | 4r? dr = An RC. 
0 0 0 


For the right-hand-side of the divergence theorem, we have u — R?f and dS = FR? sin 0 d0 d, so that 


2m pr 
pu-as = f n R* sin 0 d0 d$ = Ax R^. 
S 0 0 
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130 LECTURE 48. DIVERGENCE THEOREM (EXAMPLE 2) 


Problems for Lecture 48 
1. Test the divergence theorem using u = f/r for a sphere of radius R centered at the origin. 


2. Compute the flux integral of r = xi + yj + zk over a sphere of radius R by applying the divergence 
theorem to convert the flux integral into a volume integral. 


3. Consider the velocity field of a fluid given by 


0 A(xi t yj 4 zk) 
An (x? + y? + 22)3/2° 


u(x,y,z) 


4) Using spherical coordinates, show that 


Af 
un cau 


b) Using spherical coordinates, show that V - u = 0 provided r + 0. 
c) Using the divergence theorem, show that 


| v^ 
V 


proved that the volume V contains the origin, and is zero otherwise. You have therefore shown 
that the divergence of the velocity field is given by 


V-u=Ad(r), 


where ó(r) is the three-dimensional Dirac delta function. This velocity field is called a source 


flow. 


Solutions to the Problems 


Lecture 49 | Continuity equation 


View this lecture on YouTube 


The divergence theorem is often used to derive a continuity equation, which expresses the local con- 
servation of some physical quantity such as mass or electric charge. Here, we derive the continuity 
equation for a compressible fluid such as a gas. Let p(r, t) be the fluid density at position r and time 
t, and u(r, t) be the fluid velocity. We will assume no sources or sinks of fluid. We place a small test 
volume V in the fluid flow and consider the change in the fluid mass M inside V. 


The fluid mass M in V varys because of the mass flux through the surface S surrounding V, and 


dM 
2 = — > puds. 
dt fou m 


Now the fluid mass is given in terms of the fluid density by 


M= | pav, 
V 


and application of the divergence theorem to the surface integral results in 


a f eav =- f v- pwav. 


Taking the time derivative inside the integral on the left-hand side, and combining the two integrals 


[+v w) dV =0. 


Since this integral vanishes for any test volume placed in the fluid, the integrand itself must be zero, 


one has 


yields 


and we have derived the continuity equation 


9p 


3 ^V (eu) =0. 


For an incompressible fluid for which the density p is uniform and constant, the continuity equation 
reduces to 
V-u=0. 


A vector field with zero divergence is called incompressible or solenoidal. 
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132 LECTURE 49. CONTINUITY EQUATION 


Problems for Lecture 49 
1. Show that the continuity equation can be written as 


9 
Sy tu Vp+pV-u=0. 


2. The electric charge density (charge per unit volume) is give by p(r,t) and the volume current 
density (current per unit area) is given by J(r,t). Local conservation of charge states that the time 
rate of change of the total charge within a volume is equal to the negative of the charge flowing out of 


that volume, resulting in the equation 


5 f eena - - diras. 


From this law of charge conservation, derive the electrodynamics continuity equation. 


Solutions to the Problems 


Practice Quiz Divergence theorem 


1. The integral of u = yzi + xzj + xyk over the closed surface of a right circular cone with radius R 


and length L and base in the x-y plane is given by 
a) 0 
b) mRLVR2 +12 
c) 2n RL V R2 + 12 
d) 3TRL VR? - L2 


2. The surface integral $ r -dS over a right circular cylinder of radius R and length L is equal to 
S 
a) 0 
b) TR?L 
c) 2x R?L 
d) 3rtR?L 
3. Which velocity field is not incompressible (V - u # 0)? 
. 15. 
a) u=xyi— yj 


b) uw=(14+x)i+(1-y)j 


d) u=(x+y) + (x-y? j 


Solutions to the Practice quiz 
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134 PRACTICE QUIZ. DIVERGENCE THEOREM 


Lecture 50 | Green's theorem 


View this lecture on YouTube 


Green's theorem is a two-dimensional version of Stokes' theorem, and serves as a simpler introduction. 
Let u = u1(x,y)i + u»(x, y)j be a differentiable two-dimensional vector field defined on the x-y plane. 


Green's theorem relates an area integral over S in the plane to a line integral around C surrounding 


du» ðu u 
[ (= ai ds p(n dx + uzdy). 


this area, and is given by 


We will first prove Green's theorem for a rectangle 
with sides parallel to the axes. Let the rectangle be 
defined by a <= x <= b and c <= y <= d, as 
pictured here. The area integral is given by 


du du du d Ou 
¡(Es [on f n . 


The inner integrals can be done using the fundamental theorem of calculus, and we obtain 


d b 
f (= a) dS = / [u>(b,y) — u»(a, y)] dy «f [us (x, c) — u1(x,d)] dx = di dx + uz dy). 
Note that the line integral is done so that the bounded area is always to the left, which means coun- 
terclockwise. 

Now, given any closed smooth curve in the x-y plane enclosing an area, we can subdivide the 
area by a very fine two-dimensional rectangular grid and apply the above result to each rectangle 
in the grid. All the area integrals over the internal rectangles add. The internal rectangles share 
connecting sides over which the line integrals cancel, and the only line integrals that remain are those 
that approximate the given bounding curve. The result is Green's theorem for any area S in the plane 
bounded by a curve C. 
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136 LECTURE 50. GREEN'S THEOREM 


Problems for Lecture 50 


1. Test Green's theorem using u — —yi + xj for a square of side L lying in the first quadrant with 


vertex at the origin. 
2. Test Green's theorem using u = —yi + xj for a circle of radius R centered at the origin. 


Solutions to the Problems 


Lecture 51 Stokes’ theorem 


View this lecture on YouTube 


Green's theorem for a differentiable two-dimensional vector field u = uı(x,y)i + uz(x,y) j and a 


smooth curve C in the x-y plane surrounding an area 5 is given by 


Qu» ðu u 
i & a) as jen dx + uzdy). 


Green's theorem can be extended to three dimensions. With u = u (x,y,z) i + u2(x,y,z) j, we see that 


Qu» Qui i 
mU Y k; 


and with dS = k dS and u4 dx + u2 dy = u- dr, Green's theorem can be rewritten in the form 


] vaso p uar. 


This three-dimensional extension of Green's theorem is called Stokes' theorem. Here, 5 is a general 
three-dimensional surface bounded by a closed spatial curve C. A simple example would be a hemi- 
sphere located anywhere in space bounded by a circle. The orientation of the closed curve and the 
normal vector to the surface should follow the right-hand rule. If your fingers of your right hand point 
in the direction of the line integral, your thumb should point in the direction of the normal vector to 
the surface. 
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Problems for Lecture 51 

1. From Stokes' theorem, determine the form of Green's theorem for a curve lying in the 
a) y-z plane; 
b) z-x plane. 


2. Test Stokes” theorem using u = —yi + xj for a hemisphere of radius R with z > 0 bounded by a 


circle of radius R lying in the x-y plane with center at the origin. 


3. Consider the velocity field of a fluid given by 
| Tf-yitxj 
u(x,y) 2m ( ) ; 


a) Using polar coordinates, show that 


ró 
u(r, 0) = LT 


b) The vorticity field of the fluid is defined as w = V x u. Using polar coordinates, show that 
w = 0 provided r 0. 


c) Using Stokes theorem, show that the integral of the vorticity field over a small area in the x-y 
plane containing the origin is equal to T and therefore that the vorticity is given by w = I'ó(r), 
where ó(r) is the two-dimensional Dirac delta function in the x-y plane. This is the definition of 


a point vortex of strength T. 


Solutions to the Problems 


Practice Quiz | Stokes” theorem 


RSS 
1. Let u = —yi + xj. Compute $ u dr for the quar- > [Y 
ter circle of radius R as illustrated. Here, it is sim- 
pler to apply Stokes' theorem to compute an area " 
L > 


integral. The answer is 


a) 0 
b) inR? 
c) TR? 


d) 2nR? 


2. Letu — a 7 j. Compute the value of / (V x wu) -dS over a circle of radius R centered 


x 
x2 + y? 
at the origin in the x-y plane with normal vector k. Here, because u is singular at r = 0, it is necessary 


to apply Stokes” theorem and compute a line integral. The answer is 


3. Let u = —xlyi + xy* j. Compute $ u: dr for a unit square in the first quadrant with vertex at the 
C 


origin. Here, it is simpler to compute an area integral. The answer is 


= a 
== == 
7? WIN WIR O 


a a 
— — 


Solutions to the Practice quiz 
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140 PRACTICE QUIZ. STOKES' THEOREM 


Lecture 52 | Meaning of the divergence 
and the curl 


View this lecture on YouTube 


With u a differentiable vector field defined inside and on a smooth closed surface S enclosing a 
volume V, the divergence theorem states 


[vr wav=Pu-as. 
V 5 


We can limit this expression by shrinking the integration volume down to a point to obtain a coordinate- 
free representation of the divergence as 


v0 


1 
V -u = lim y Qu ds. 


Picture V as the volume of a small sphere with surface S and u as the velocity field of some fluid of 
constant density. Then if the flow of fluid into the sphere is equal to the flow of fluid out of the sphere, 
the surface integral will be zero and V - u = 0. However, if more fluid flows out of the sphere than in, 
then V -u > 0 and if more fluid flows in than out, V -u < 0. Positive divergence indicates a source 
of fluid and negative divergence indicates a sink of fluid. 

Now consider a surface S bounded by a curve C. Stokes” theorem states that 


[ vaso fuam 


We can limit this expression by shrinking the integration surface down to a point. With n a unit 
normal vector to the surface, with direction given by the right-hand rule, we obtain 


(Od) esi t u-dr. 
s>0 S É 


Picture S as the area of a small disk bounded by a circle C and again picture u as the velocity field 
of a fluid. The line integral of u - dr around the circle C is called the flow's circulation and measures 
the swirl of the fluid around the center of the circle. The vector field w = V x u is called the vorticity 
of the fluid. The vorticity is most decidedly nonzero in a wirling (say, turbulent) fluid, composed of 
eddies of all different sizes. 
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142 LECTURE 52. MEANING OF THE DIVERGENCE AND THE CURL 


Problems for Lecture 52 
1. The incompressible Navier-Stokes equation governing fluid flow is given by 


ou = 1 2 
a ev oe u, 


with V-u = 0. Here, p and v are fluid density and viscosity. 


4) By taking the divergence of the Navier-Stokes equation, derive the following equation for the 


pressure in terms of the velocity field: 


b) By taking the curl of the Navier-Stokes equation, and defining the vorticity as w = V x u, derive 


the vorticity equation 


9 
P + (u- V)w = (w: V)u 4 vV?w. 
You can use all the vector identities presented in these lecture notes, but you will need to prove 


that í 
u x(V xu)= zV (u-u) - (u: V)u. 


Solutions to the Problems 


Lecture 53 | Maxwell's equations 


View this lecture on YouTube 


Maxwell's equations in SI units and in integral form are given by 


$ E-dS = dem, (Gauss's law for electric fields) 

S 0 

$ B.dS —0, (Gauss's law for magnetic fields) 
S 

$ E - dr = E i. B-dSs, (Faraday’s law) 
E dt Js 

$ B -dr = pọ (I. + 0 f E: is) A (Ampere-Maxwell law) 
C S 


where E and B are the electric and magnetic fields, genc and Tenc are the charge or current enclosed 
by the bounding surface or curve, and £o and uy are dimensional constants called the permittivity and 
permeability of free space. 

The transformation from integral to differential form is a straightforward application of both the 
divergence and Stokes” theorem. The charge q and the current I are related to the charge density p 


and the current density J by 
a= | pav, r= f 3-48. 
V S 


We apply the divergence theorem to the surface integrals and Stokes' theorem to the line integrals, 
replace Jenc and lenc by integrals, and combine results into single integrals to obtain 


Io 8-2) dV =0, | (vs 252) as- o 
V £0 E dt 
JE 
], v mav-o [(v«5 n (rro )) as =0 
- A ot 


Since the integration volumes and surfaces are of arbitrary size and shape, the integrands must vanish 


and we obtain the aesthetically appealing differential forms for Maxwell’s equations given by 


V.E- P, Vere 
£0 ot 
OE 
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144 LECTURE 53. MAXWELL'S EQUATIONS 


Problems for Lecture 53 


1. Using Gauss's law for the electric field given by 
$ E-ds = dee, 
S E0 


determine the electric field of a point charge q at the origin. Assume the electric field is spherically 


symmetric. 
2. Using Ampère’s law given by 


B dr = Holenc, 
Cc 


determine the magnetic field of a current carrying infinite wire placed on the z-axis. Assume the 


magnetic field has cylindrical symmetry. 


Solutions to the Problems 


Appendix 
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Appendix A | Matrix addition and 


multiplication 


View this lecture on YouTube 


Two-by-two matrices A and B, with two rows and two columns, can be written as 


b 
pe [n 
boy 


b» 
boo 


The first row of matrix A has elements a11 and 412; the second row has elements a2; and a». The 


first column has elements 41, and 421; the second column has elements a12 and a». Matrices can be 


multiplied by scalars and added. This is done element-by-element as follows: 


k k 
ka= |u "2| A+B= 
kazı kan 


ay bg an 
ay, b» an 


tb) 
tb» 


Matrices can also be multiplied. Matrix multiplication does not commute, and two matrices can be 


multiplied only if the number of columns of the matrix on the left equals the number of rows of the 


matrix on the right. One multiplies matrices by going across the rows of the first matrix and down the 


columns of the second matrix. The two-by-two example is given by 


an ay (bu ba) 
an an) Aba ba 


ayıbıı 4 


+ a12b21 


a21b11 4 


+ a22b21 


d11012 + 


azıbı2 + 


+ 1252» 
+ a22b22 


Making use of the definition of matrix multiplication, a system of linear equations can be written 


in matrix form. For instance, a general system with two equations and two unknowns is given by 


81131 + 419%) = by, 


and the matrix form of this equation is given by 


In short, this matrix equation is commonly writt 


A21Xı + d22X2 = bo; 


en as 


Ax — b. 
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148 APPENDIX A. MATRIX ADDITION AND MULTIPLICATION 


Appendix B | Matrix determinants and 


inverses 


View this lecture on YouTube 


We denote the inverse of an n-by-n matrix A as A^, where 


AAT! = AİA =], 


and where I is the n-by-n identity matrix satisfying IA = AI = A. In particular, if A is an invertible 
matrix, then the unique solution to the matrix equation Ax = b is given by x = Ab. 

It can be shown that a matrix A is invertible if and only if its determinant is not zero. Here, we only 
need two-by-two and three-by-three determinants. The two-by-two determinant, using the vertical bar 
notation, is given by 


a1 412 
= 411422 — 012051; 


421 422 


that is, multiply the diagonal elements and subtract the product of the off-diagonal elements. 
The three-by-three determinant is given in terms of two-by-two determinants as 


aa An 413 
421 422 423 = 411 


The rule here is to go across the first row of the matrix, multiplying each element in the row by the 
determinant of the matrix obtained by crossing out that element’s row and column, and adding the 
results with alternating signs. 

We will need to invert two-by-two and three-by-three matrices, but this will mainly be simple 
because our matrices will be orthogonal. The rows (or columns) of an orthogonal matrix, considered 


as components of a vector, are orthonormal. For example, the following two matrices are orthogonal 


matrices: . : i 
SEE . sinÜcos$ sin@sing cos 0 
cosÜ sinó . . 
. j cosücos$ cos@sing  —sin0 
—sin@ cos® . 
— sin $ cos $ 0 
For the first matrix, the row vectors f = cos ĝi + sin0j and 0 = — sin ĝi + cos0j have unit length and 


are orthogonal, and the same can be said for the rows of the second matrix. 
The inverse of an orthogonal matrix is simply given by its transpose, obtained by interchanging 


the matrices rows and columns. For example, 
-1 
cosÜ sind | [cos0 —siné 
—sin@ cos0 sin 9 cos]. 


For more general two-by-two matrices, the inverse can be found from 


i.e., switch the diagonal elements, negate the off-diagonal elements, and divide by the determinant. 
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150 APPENDIX B. MATRIX DETERMINANTS AND INVERSES 


Appendix C | Problem and practice quiz 


solutions 


Solutions to the Problems for Lecture 1 


1. We show the associative law graphically: 


(A+B)+C=A+4+(B+C) 


2. Draw a triangle with sides composed of the vectors A, B, and C, with C = A+ B. Then draw the 
vector X pointing from the midpoint of C to the midpoint of B. 
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From the figure, we see that 


1 1 
pC TX - Ac 5B. 


Using C — A 4 B, this equation becomes 


1 
(A-B)*X-A-3B, 


NI =e 


and solving for X yields X = 1A. Therefore X is parallel to A and one-half its length. 


Solutions to the Problems for Lecture 2 


a) The unit vector that points from mı to mz is given by 


mor . (x2—x1)i+ (y2 —y1)j + (z2 — zı)k 
[ro — "1| V (x2 = 1)? + (ya — y1)? + (22 — 21)? 


b) The force acting on mı with position vector rı due to the mass mz with position vector rz is 


written as 


(xa — x1)d + (y2 - y1)j + (z2 — zi)k 
[(xo — x1)? + (yo — y1)? + (22 — 21)2]9/2" 


T5) —T1 
In — rıl? 


F = Gmıma = Gmim» 


Solutions to the Problems for Lecture 3 


a) A- B = AıBı + A5B5 + A3B3 = B141 + B242+8B343 = B- A; 


b) A-(B+C) = A1(B4 + C1) + A»(B5 + C2) + Az(B3 + C3) = A,B, + A4C1 + A2B + A2C5 4 
A3B3 + A3C3 = (A1B1 + A2B2 A3B3) | (A1C4 + AoCo + A3C3) — A- B- A: C; 


c) A. (kB) = Aj (kB1) + A2(kB2) + A3(kB3) = (kA1)B1 + (kA2)Ba + (kA3) Ba = k(A1B1) +k(A2B2) + 
k(A3B3) = (kA)- B — k(A- B) 


2. The dot product of a unit vector with itself is one, and the dot product of a unit vector with one 


perpendicular to itself is zero. That is, 


ivi=j.j=k k=l; d joi kk-jk-0; j-i=k-i=k-j=0. 


3. Consider the triangle composed of three vectors pictured below. 
With C — A — B, we have 


IC? 2 C-C 2(A- B). (A- BJ) 2 A A- B. B-2A. B - |AP +/|B/? —2|A||B|cos0, 
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where 0 is the angle between vectors A and B. In the usual notation, if A, B and C are the lengths of 


the sides of a triangle, and 0 is the angle opposite side C, then 


C? = A? + B? — 2ABcos0. 


Solutions to the Problems for Lecture 4 


1. 
a) 
i j k j k 
AxB= A4 A» As | = — | Bı Bo B3 =-BxA. 
Bı B» Bs A, An As 
b) 
i j k 
Ax (B + C) = Ay A2 A3 
By +C, BotCon B3+C3 
i j k i j k 
=|A, An Al -|A; An 43 - Ax B Ax C. 
Bı B; B| |C & C3 
c) 
i j k i j k i j k 
Ax (kB) = Ay A» Az =k Ay A) Az = kAı kA5 kA3 
kB, kB kB; Bi Bo Bs Bi de B 


— k(A x B) = (kA) x B. 


2. The cross product of a unit vector with itself is equal to the zero vector, the cross product of a unit 


vector with another (keeping the order cyclical in 2, j, k) is equal to the third unit vector, and reversing 
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the order of multiplication changes the sign. That is, 
ixi=0, jxj=0, kx k= 0; 


ixj=k, jxk=i, kxi=j; 


kxj=-i, jxi=-k ixk=-j. 


3. One such example is 


ix(ixk)=-ixj=-k, 
(ixi)xk=0xk=0. 


Solutions to the Practice quiz: Vectors 


1. c. As an example, ix (ix j) + (i x i) xj. 


2. b. 
i j k 
(Ax B):.j— am a a3|:j = 43b, — abs. 
7 
3. d. 


ix(jxk)=ixi=0, (ixj)xk=kxk=0, (ixi)xj=0xj=0, ix(ixj)=ixk=-j. 


Solutions to the Problems for Lecture 5 


1. We first compute the displacement vector between (1,1,1) and (2,3,2): 


u = (2-—1)t+ (3 —1)j + (2—1)k = i +2j +k. 
Choosing a point on the line to be (1,1,1), the parametric equation for the line is given by 


r= ro+ut = (i+j+k)+(i+2j+k)t= (1+ t)i +(1+2t)j +(1+t)k. 


The line crosses the x = 0 and z = 0 planes when t = —1 at the intersection point (0, —1,0), and 
crosses the y = 0 plane when t = —1/2 at the intersection point (1/2,0,1/2). 
Solutions to the Problems for Lecture 6 


1. We find two vectors parallel to the plane defined by the three points, (—1, —1, —1), (1,1,1), and 
(1, 1,0): 


sı = (141)i+(141)9 + (141)k = 26427 42k, — 8; = (1-1) (-1—1)j  (0—1)k = —2j — k. 
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We can divide s, by 2 to construct a normal vector from 


i i j k 
N=;5sx2=|1 1 1|=i+j-2k 
02 zd 


The equation for the plane can be found from N - (r — r2) = 0, or IN. r = NN : ro, or 
(i--j —2k)- (xi yj zk) =(4+j—-2k)-(i+g7 +k), or xty—2z=0. 


The intersection of this plane with the z — 0 plane forms the line given by y — —x. 


Solutions to the Practice quiz: Analytic geometry 


1. d. Write the parametric equation as r = ro + ut. Using the point (0, 1,1), we take ro = j + k and 
from both points (0,1,1) and (1,0, —1), we have u = (1-0): + (0—1)j + (-1—1)k = i — j — 2k. 
Therefore r = j +k + (i— j — 2k)t = ti + (1 — t)j + (1 — 2t)k. 


2. a. The line is parameterized as r = ti + (1 — t)j + (1— 2t)k. The intersection with the z = 0 plane 
1 1 11 
occurs when t = 1/2 so that r = 5! + 2 j. The intersection point is therefore (5; > 0). 


3. d. We first find the parametric equation for the plane. From the points (1,1,1), (1,1,2) and (2,1,1), 
we construct the two displacement vectors 


The normal vector to the plane can be found from 
N=3s xs. =kx(i-k)=kxi-kxk=j. 


Therefore, the parametric equation for the plane, given by N - (r — r1) — 0, is determined to be 


or y = 1. This plane is parallel to the x-z plane and when z = 0 is simply the line y = 1 for all values 
of x. Note now that we could have guessed this result because all three points defining the plane are 
located at y = 1. 


Solutions to the Problems for Lecture 7 


a) If ijk is a cyclic permutation of (1,2,3), then ejj = Ejki = €xij = 1. If ijk is an anticyclic per- 


mutation of (1,2,3), then ei; €jki €kij —1. And if any two indices are equal, then 


Eijk = Ei = €kij = 0. The use is that we can cyclically permute the indices of the Levi-Civita 
tensor without changing its value. 
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b) If ijk is a cyclic permutation of (1,2,3), then ej, = 1 and ej; = €kji = eg; = —1. If ijk is an 
anticyclic permutation of (1,2,3), then ej = —1 and ejy = €xji = ei; = 1. And if any two 
indices are equal, then €jjx = €jik = €ikj = 0. The use is that we can swap any two indices of the 
Levi-Civita symbol if we change its sign. 


2. We have 
E1jk AjBk = €12342B3 + €13243B2 = A2B3 — A3B = [A x B] 
E2jk AjBy = €231A3By + €21341B3 = AzB1 — A1B3 = [A x B], 
€ajkAjBy = €31241B2 + €32142B1 = A1B? — A2B1 = [A x Bla. 
3. 


a) Now, 6;;Aj = 91141 + 0242 + 0343. The only nonzero term has the index of A equal to i, 
therefore ¡Aj = Aj. 


b) Now, ojxdkj = 0101; + i252; + 01363). If 1 + j, then every term in the sum is zero. If i = j, then only 
one term is nonzero and equal to one. Therefore, ö;köx; = ójj. This result could also be viewed as 
an application of Part (a). 


4. We make use of the identities ô; = 3 and ójjój, = ójj. For the Kronecker delta, the order of the 
indices doesn’t matter. We also use 


€iik€Imn = Óil(ÉjmÓkn — OjnOkm) — Sim (den — Ojndk1) + Fin (9j19km — OjmÓxi)- 


Ciik€imn = Oii(SjmOkn — ÓjnÓkm) — Sim (OjiOkn — Oindki) + Fin (Sj Okm — Ojm Ki) 
= 3(dimÓkn — 9jn%km) — (OjmÓkn — ÓjnÓkm) + (9jnkm — ÓjmÓka) 
= dimOkn = ÜinÓkm- 


b) We use the result of a) and find 


€ijk€ijn = Öjjökn m din Oki = 30kn — Ôkn = 20kn- 


Solutions to the Problems for Lecture 8 


1. We sometimes need parentheses because the vector product is not associative and expressions can 


be evaluated in more than one way with different results. We can write without any ambiguity, 


(scalar triple product) A-BxC=B:-CxA=C-AxB, 
(vector triple product) Ax(BxC)=A:CB-A:BC 
(scalar quadruple product) AxB-CxD=A:CB-D-A-DB-C 
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For the vector quadruple product, we can write without ambiguity 
(Ax B) x (Cx D) = (Ax B). DC-(Ax B). CD 
or 
(AxB)x(CxD)=(AxB-D)C-(AxB-C)D. 


The absence of any parentheses on the right-hand side would result in an ambiguity, however, since 
in general, 
(AxB):D)C%A x ((B-D)C). 


The vector on the left-hand side of this expression is parallel to C while the vector on the right-hand 
side is perpendicular to C. The two expressions can therefore be only equal if C is the zero vector. 


Solutions to the Problems for Lecture 9 


1. Consider A- (B x C). If A = B, then A: (Bx C) = A. (Ax C) =0 since A x C is orthogonal 
to A. A similar results holds for A = C. If B = C, then A- (Bx C) = A: (B x B) = 0 since 
BxB=0. 


2. We use the fact that the value of the scalar triple product is unchanged under a cyclic permutation 
of the three vectors, and that the dot product is commutative. We have 


A-BxC—C.Ax B—Ax B.C. 


3. Consider the scalar triple product 
ej: (ej x ey). 


We have already proved that if any two indices are equal, then the scalar triple product is zero. 
Furthermore, the scalar triple product is unchanged under a cyclic permutation of the three vectors so 
that 


€i: (e x e3) = e: (ez x e1) = ea: (ez x e2), 


and 
e1: (e2 x es) =i- (jxk)=i- i=1. 
Finally, 
e3° (ez x e1) =e2: (e x e3) =e1: (ez x e2), 
and 


e3: (e2 x e) =k- (j x i) 2 -k-k = —1. 


We have computed all possible cases, and have thus proved 
ej: (ej x ex) = Eijk- 


A slicker proof would denote the mth component of e; by ejm and use ejm = djm, and so on for the 
other unit vectors and components. Then 


ej: (ej x ey) = Eil€Imn€jm€kn = Elmndi9jmÚdkn = Ejjk- 
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Solutions to the Problems for Lecture 10 


1. We prove the Jacobi identity: 


Ax(BxC)-Bx(Cx A) -Cx(Ax B) 

= [((A- C)B - (A - B)C] + ((B- A)C - (B- C)A] + [(C- B)A — (C- A)B] (vector triple product) 
—[(A-C)B - (C AJB] -|(B- A)C - (A: B)OC] x ((C- BJA- (B: C)A] (rearranging terms) 
=0. 


2. We want to prove that the scalar quadruple product satisfies 


(Ax B): (C x D) - (A. C(B- D) - (4. D(B.C). 


We have 
(A x B)- (C x D) = [A x B];[C x DJ; (X-Y 2 XY) 
= Ej AjBr€iimC¡Dm (X x Y |; = eig X; Y) 
= eijk€itm Aj BkCi Dm 
= (Ój1Ókm — Ójm0x1) AjBxCiDm (Eijk€ilm = Ójlóxm — Ójmóki) 
= AjC;B,D, — AjDjByC (S/C, = Cj, etc.) 
=(A-C)(B-D)—(A-D)(B-C). (A/C; = A- C, etc) 


3. We can prove Lagrange's identity using the scalar quadruple product identity (A x B): (Cx D) = 
(A: C)(B- D) —(A-D)(B-C). We have 


|Ax B (Ax B). (Ax B) - (A A(B. B) - (A B(B- A) - |AP|B? - (A- By. 
An alternative proof uses 


|A x Bl? = |AP|Bp? sin? 0 = | AP|B[?(1— cos? 0) = | AP|BP? — | AP| B cos” 6 = | AP|BP? — (A- By. 


4. We want to prove that the vector quadruple product satisfies 
(Ax B) x (Cx D)  ((Ax B): D)C — ((Ax B): C)D. 
We will make use of the vector triple product identity given by 
Ax(BxC)=(A-C)B-(A-B)C. 
Let X = A x B. Then using the vector triple product identity, we have 


(Ax B)x (Cx D) =X x (Cx D) 
=(X:D)C-(X:C)D 
— ((Ax B): D)C — ((Ax B): C)D. 
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Solutions to the Practice quiz: Vector algebra 


1. c. The relevant formula from the lecture is €;¡x€¡lm = ójjógg — 9jmOx1- To directly apply this formula, 
we permute the indices of the Levi-Civita symbols without changing their cyclic order: 


€ijk€ljm = €jki€jml = OkmOil — Óx1Ójm- 


2. d. The other expressions can be shown to be false using A x B — —B x A and, in general, 
Ax(BxC)*(Ax B)xC. 


3. c. Use the facts that A x B is orthogonal to both A and B, A- B is zero if A and B are orthogonal, 
and A x B is zero if A and B are parallel. 


Solutions to the Problems for Lecture 11 


1. Scalar fields: electrostatic potential, gravitational potential, temperature, humidity, concentration, 
density, pressure, wavefunction of quantum mechanics. 
Vector fields: electric and magnetic fields, magnetic vector potential, velocity, force fields such as gravity. 


Solutions to the Problems for Lecture 12 


1. Using the chain rule, 


of —2nx of —2ny of —2nz 


Ox (x2 + y? +22)n+1" dy (x? + y? + 22)n+1’ oz (x2 + y2 + 22)n+1" 


2. Define 
flt+e,x +0) = gle, ô). 


Then the first-order Taylor series expansion of g is given by 
g(€,5) = g(0,0) + eg+(0,0) + dgx(0,0), 
which in terms of f becomes 
flit+ex+0)=f[t,x)+efilt,x) +ôfx(t, x). 
Applying this expansion to f(t + «At, x + BAtf (t, x)), we have to first-order in At, 


f(t -- &At,x + BAtf(t,x)) = f(t, x) +aAtfı(t,x) + BAtf(tx) (t, x). 
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Solutions to the Problems for Lecture 13 


1. The formulas derived in the text are 


| ox Dy X xwiXxi 
ny — (Zx)? 


_ ay xii — xu) y) 
nye -— (Lx)? 


Bo , Bi 


, 


where the sum is from i = 1 to 3. Here, x1 1, x? = 2, x3 = 3, and y; = 1, y2 = 3 and y3 = 2. We 


_ (14)(6) — (13)(6) _ _ 93) = (6)(6) _ 
O A P- Gag-er A 


The best fit line is therefore y = 1 + x/2. The graph of the data and the line are shown below. 


have 


T T T 
3 ° =i 
c 2 ° " 
1 e J 
| i | 
1 2 3 
XL 


Solutions to the Problems for Lecture 14 
1. 


a) With f(x,y) = e™, x = rcos0, and y = rsin 0, application of the chain rule results in 


of  ofox  ofoy 
dr oxor  Oyor 


= ye"! cos 0 + xe™ sin 


1? cos sin 1? cos 0 sin 0 


= r sin Ó cos de +rsin 0 cos 0e 


: 2 i 
= 2r sin 0 cos be" <os0sin®, 
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and 
af  afax ,Of9y 
90  Ox00 ay ae 
= ye™ (—r sin 0) + xe™ (r cos 0) 


2 r? 


a . E . 
= -r sin? 0e cossint | „2 cos? ger cos 0 sin 0 


; 2 i 
= r?° (cos? 0 — sin? gg ce6sne, 


b) Substituting for x and y, we have f = e” cosósinó. Then 


i = 2r cos sin be" «95? Sind, 
or 
of = r? (cos? 0— sin? 8)e" cos @ sin 6 


Solutions to the Problems for Lecture 16 
1. Suppose ax + by + cz = 0. We have the relations 


_ —by -= cz TAX — cz _ —ax — by 
xm I A acr cS EC 


x 


The partial derivatives are 


and the triple product is 


ES 


2. Suppose ax + by + cz + dt = 0. We have the relations 


ELS — cz — dt —ax — cz — dt „_ Zr by dt 


. —ax — by — cz 
a d b i c : d f 


x 


The partial derivatives are 


and the quadruple product is 
Ox dy dz dt _ b C d ay 
dy dz ot dx ( =) ( 3 ( =) ( a) n 


Apparently an odd number of products yields —1 and an even number of products yields 4-1. 


Solutions to the Practice quiz: Partial derivatives 


1. d. The partial derivative with respect to x is given by 


Of _ x 
Ox (x2 + y2 + 22)3/2’ 
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and the mixed second partial derivative is then given by 


of u 3xy 
oxoy (x? +y2 +22)5/2 


2. a. From the data points (0,1), (1,3), (2,3) and (3,4), we compute 


Yo=6 33314 )y 1L Y xy =21 


Then using 


| Ou Ey Y 58 
Axa = De 


miu -—QO xy) 


LAS nz 


Bo 


we have 


_ (14)(11) — (21)(6) 154—126 28 7 


Po = ps) — (6)? 56—36 20 5 
(4)(21)—(6)(11) 84-66 18 9 


Pi = aaa) (6) ~ 56-36 20 10 


The least-squares line is therefore given by y = 7/5 + 9x/10. 


3. d. Let f = f(x, y) with x = rcos0 and y = rsin@. Then application of the chain rule results in 


df Ofox , af dy 


900 09x00 dy 00 
. „of i of 
= reino | rcos = 
__ NF 


Solutions to the Problems for Lecture 17 


a) Let p(x,y,z) = x? + y? + 2?. The gradient is given by 
Vo = V (X +4? +27) = 2xi + 2yj +2zk. 
In terms of the position vector, we have 
V(r?) = 2r. 


1 
VA 


b) Let p(x,y,z) . The gradient is given by 


1 
Vo=V 
? (+=) 
x z 


; y ; 
k. 
(x? 4 y2 + 22)872° HFI 2 y y2 4 22/372 
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In terms of the position vector, we have 


Solutions to the Problems for Lecture 18 


a) With F = xyi + yzj + zxk, we have 


Solutions to the Problems for Lecture 19 


a) With F = xyi + yzj + zxk, we have 


i j k 
VxF=l0/0x ð/ðy 90/0z| = —yi —zj — xk. 
xy yz ZX 


b) With F = yzi + xzj + xyk, we have 


i j k 
Vx F=|0/dx 9/dy 0/02 = (x—-x)it (y—y)3 + (z-z)k = 0. 
yz xz — xy 


2. With u = u1(x, y)? + ua(x,y)j, we have 


i j k 
— = EE m Qu» ðu = 
w=Vxu=! 0/9 0/0y  0/O0z| =i +0j 4 (5: GE w3k. 


u(x,y) ulxy) 0 


Therefore, w3 = du) /d0x — du1 / dy. 
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Solutions to the Problems for Lecture 20 


1. We have 


(1 1 a» 1 90) 1 
r ax? | V+ +Z) ol ytz) 32 | rgrz]' 


We can compute the derivatives with respect to x and use symmetry to find the other two terms. We 


have 
b 1 i =x . 
ox [x2 + y2 +22 (x? + y? + 22)3/2’ 
and 
dx \ (x2 + y2 + 22)3/2 (x? 42 +22)? 
1 | 3x? 
(x2 + y2 + 22)3/2 | (x2 + y2 + 22)5/2° 


It is easy to guess the derivative with respect to y and z, and we have 


v (2) = 3 | 3(x? + y? +2?) 
T (x2 + y2 + 22)3/2 * (124 y2 + 22)5/2 
3 3 


(x2 + y2 + 22)372 (x2 + y2 + 22)3/2 = 0, 


a result only valid for r 4 0. 


Solutions to the Practice quiz: The del operator 


1. d. We have 
1 1 
v(a) =v (rr) 
7 —2x i —2y " —2z 
CER FFA prey 

a 2r 

x 
2. b. We use 

Y.F-V. xi+yj+zk 
va? + y? + 22 ` 

Now, 


ð x 5 VE y? z2 x(x? i y? | ae 
ox V Y2 Fy +22 2+y+22 


1 x? 
r y 


and similarly for the partial derivatives with respect to y and z. Adding all three partial derivatives 


results in F 5 5 
v.p-? a _ 3 1.2 
r T T r r 
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3. 
b. We have 
i j k 
Vxr=|9/dx 0/dy 0/0z| =Q. 
x y zZ 


Solutions to the Problems for Lecture 22 


a) To prove V - (fu) =u- Vf - fV: u, we compute 


V (fu) = 2 (fui) (divergence in component notation) 
1 
_ of Qui - 
= dc^ itf df (product rule for the derivative) 
=u Vf+fV-u. (back to vector notation) 


b) To prove V x (V x u) = V(V - u) — V?u, we compute for the ith component: 


d 9 
[V x(V x u)); = Eijk (eun =) (curl in component notation) 
ox; Ox] 
um " 
= Éiik€kIn Ix axi (Ekim doesn't depend on xj) 
— Ene um (Eik = Ekij) 
= EkijEklm Oxjoxj ijk = ©kij 
um 
= (óiójm — LETT (Ekij€kim = Sirdjm — Simdjr) 
Pu; du; " OF th _ dum ai 
Oxjüx;  OxjOx; i Oxj0x] rar 
= [V(V - u)]; - [V?u];. (back to vector notation) 
Therefore, V x (V x u) = V (V -u) — V?u. 
2. 
a) With dr/dt = u(t,r(t)), the component equations are given by 
dy dx dx 
a = u1 (t; x1, X2, Xa), = = Up (t; x1, X2, X3), D = ua(t; x1, x2, X3). 
dt dt dt 


b) Using the chain rule, we can compute the second derivative of x, as 


dx; ou ðu dxı i ou dx» j du; dxa 
dt? ot od om dt 9x3 dt 
= Q4 i du; i du; i du; 

a aa Caeo E 
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Similarly for x? and xs, we have 


dx, = Qu» Qu» T Qu» ón du) d’xz E Qus er Qus dy Qua - Qua 
e ni Can Pan da dB oy, Cun "5n 
c) Using the operator 
u V=u Fu Le u 
Dx i ?9x; | 33x3” 


the three components can be combined into the vector expression 


dr ðu 


This expression is called the material acceleration, and is found in the Navier-Stokes equation of 


fluid mechanics. 


Solutions to the Problems for Lecture 23 


1. Start with Maxwell's equations: 


V. E-0, V.B=0, Vx E=-%2, V x B = moon. 


Take the curl of the fourth Maxwell’s equation, and commute the time and space derivatives to obtain 
0 
Vx(VxB)= Ho€o a, (V x E). 
Apply the curl of the curl identity to obtain 
2 ð 
V(V-B)-V’B= Ho€o 5, (V x E). 


Apply the second Maxwell's equation to the left-hand-side, and the third Maxwell's equation to the 


right-hand-side. Rearranging terms, we obtain the three-dimensional wave equation given by 


_ 202 
ga ^"^" VB 


where c =1/,/H0€0- 


Solutions to the Practice quiz: Vector calculus algebra 


1. a. We make use of the vector identity 
V(u:v) 2 (u- Vv - (v. V)u+ux (V xv)+v x(V xu). 


Setting v — u, we have 
V(u-u) —2(u: V)u -2u x (V x u). 


Therefore, 
5V (uu) =ux(Vxu)+(u-V)u. 
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2. c. The curl of a gradient (a. and d.) and the divergence of a curl (b.) are zero. The divergence of a 


gradient (c) is the Laplacian and is not always zero. 


3. b. With E(r,t) = sin(z — ct)i, we have V - E = 0, and 
i j k 
VxE=|  0/0x 9/dy 9/9z| = cos (z — ct)3. 
sin(z—ct) 0 0 


oB 
Maxwell's equation V x E — — E then results in 


oB 
ot 


= — cos (z — ct)3, 
which can be integrated (setting the constant to zero) to obtain 


1 
B= z Sin (z — ct) J. 


Solutions to the Problems for Lecture 24 
1. In general, the mass of a solid with mass density p = p(x, y, z) is given by 


M= [|] p(,1,2)axayaz 


To determine the mass of the cube, we place our coordinate system so that one corner of the cube is at 
the origin and the adjacent corners are on the positive x, y and z axes. We assume that the density of 


the cube is only a function of z, with 


plz) =p ++ = (02—p1). 


The mass of the cube is then given by 


M= ITTI "ES (p2 — p1) | ético f ax f av f | pide = (p2—pr)| d z 


=r? [pra + kar = [p = 2 e ls +3502- m) - 5L (p t p2). 
Solutions to the Problems for Lecture 25 


The figure illustrates the integral over x first and y second. With this order, the integral over the 


parallelogram is given by 
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i 7 
z= y/3j / 
| / J 
/ 
/ 
0.6 | " 
> / / 
/ 
04+ f / 
/ / 
/ 
02 | / / 
i /e=1+y/3 
0 1 1 1 L L 
0 0.2 0.4 0.6 0.8 1 1.2 
x 
1 pl+y/3 1 ,3,, ¡1=1+y/3 
f / ” x?ydxdy =| nif dy 
0 Jy/3 o 9 xy 


1 fi jx — yq Xs 

p = E Misa d 
(e) - (9) ) 
1 1 


Solutions to the Practice quiz: Multidimensional integration 


1. b. To find the volume, we integrate z = xy over its base. We have 
1 pl 1 1 1 2 24 
/ / xydxdy= | xdx | y dy = / xdx| = (3) Lc. 
o Jo 0 0 0 2 4 


2. b. To determine the mass of the cube, we place our coordinate system so that one corner of the 
cube is at the origin and the adjacent corners are on the positive x, y and z axes. We assume that the 
density of the cube is only a function of z, with 


p(z) = (142) g/cm*. 


The mass of the cube in grams is then given by 


1 pl pl 1 1 1 a4 
mM=/ 1 n (1 +2)dedy dz = | ax | dy | (1+2) dz = (z+ 527)|, = 15g. 
o Jo Jo 0 0 0 2 


3. d. We draw a picture of the triangle and illustrate the chosen direction of integration. 


0.5[ 


Integrating first along x and then along y, the volume is given by 


1 p2-y 14 à 1 1 j j 
"M xydxdy= f z“ yay => | y|(2-y) - | dy 
0 Jy o 2 2 Jo 


; E di d 
= ze. = = 
=2 f y )dy 2 (5 5) 5 


2-y 
y 


Solutions to the Problems for Lecture 26 


1. The matrix form for the relationship between f, Ó and i, j is given by 


PY cosÜ sinó i 
6} \-sin® cos jl 


Inverting the two-by-two matrix, we have 
i| [cos0 —sin0| [f 
j — Asin cos 0 6) 


i = cos 0f — sin 06, j = sin Or + cos 06. 


Therefore, 


a) We use 
x —rcos6, y = rsin0. 
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With f = f(x(r,0), y(r,0)), we have using the chain rule, 


Qf fox | af ay Fo. „ðf 
= = cos0 2. + sing, 
or oxor  oyor gc em oy 
Qf afax  afay | af af 
a. SC 


b) We can write the result of Part (a) in matrix form as 


of /or u cos 0 sin 0 of /àx 
9f/30] I—rsing rcos0) \af/ay) 


Inverting the two-by-two matrix results in 


of /ox _1 rcos@ —sin0 of /or 
9f/9y]  rVXrsinü cose } V0f/90] ' 


Or 
of Of sindof of . „ðf , cos0 0f 
= = 0 | . 
p V Lu au r 
The Cartesian partial derivatives in polar form are therefore 
u o9 sine ð ð ing? , cos ð 
ax arr 90 Oy "Or r 0 


3. We have 
rf —rcos0i-Frsin0j = xi +yj, 


and 
rô = -rsindi+rcos0j = —yi x3. 


Solutions to the Problems for Lecture 27 


1. We have 1 
u= 7 (kf + k;6) = uy? +u90, 
so that 
"PE kı "m k2 
T = r , g — r * 
The divergence is given by 
1 0 1 dug 
Mg ra "+ 00 
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and the curl is given by 


“CEL (2)) 23, (2) 


It is important to emphasis that these results are only valid for r + 0. The vector field is singular when 
r — and is not differentiable in the usual sense. 


Solutions to the Problems for Lecture 28 


1. With u = u(r), we solve 
G 


vp' 

with boundary condition u(R) — 0. Writing the Laplacian in polar coordinates, we have 
ld „au | G 
rdr N dr vp 


We multiply by r and integrate from 0 to r: 


y T 
l A CF) ir [ rdr, 
o dr \ dr ve Jo 


V^ = — 


or 


or 


Applying the no-slip boundary condition at r = R, we obtain 


un = SE (a-(2)’). 


The maximum velocity occurs in the center of the pipe and is given by 


so we can write the velocity field as 
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Solutions to the Problems for Lecture 29 


1. We have 


I=rxp=rx (mr) = mr? x (ré +706) = mr Òf x ô). 


Now, # and Ó are perpendicular unit vectors so that |f x 0| = 1, and 


| = mr? jle x 0| = mr?|6]. 


2. Using the Einstein summation convention, we have 


d (1 2 1 d 
dt (nio ) E 5 4, 00) 


Solutions to the Problems for Lecture 31 


1. The mass density of the disk is given by 


a(r) = po + (pı — po) (r/R). 


Integrating the mass density in polar coordinates to find the total mass of the disk, we have 


27 R 
M= f [ (po-+(p1—po)(r/R)] raras 


(pı — or] m 
r=0 


2 3R 


1 
¿TR (po +2p1), 


2. It is simplest to do this integral by transforming to polar coordinates. With x? + y? = r? and 
dx dy = r dr d0, we have 


on co 5 27t co 2 00 2 
p =/ / e" rdrdó =; do | re" dr = zu | re" dr. 
0 0 0 0 0 


Let u =r? and du = 2r dr. Then the integral transforms to 


Therefore, 
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Solutions to the Practice quiz: Polar coordinates 
1. d. Using 6 = — sin di + cos 05, we have 


rÔ = —r sin bi --rcos0j = —yi + xj. 


2. b. With f = cos 01 + sin 0j and 6 = — sin bi + cos 05, we have 


dé E a 
a — cos 0 — sin 0j = —f. 


3. The mass density of the disk is given in polar coordinates by 
c = c(r) = (10 — 9r) g/cm?. 


The mass is found by integrating in polar coordinates using dx dy = r dr d0. Calculating in grams, we 


27 yl 
M= / / (10 — 9r)r dr d0 
0 0 


270 1 
= do | (10 — 9r)r dr 
0 0 


1 


have 


= 2n(5r? -3r°)| = 4 = 12.57g. 


0 


Solutions to the Problems for Lecture 32 


1. We have 
V ex +05 +22 
= (cos pp — sin $$) (cos Z^ ms 5) + (sin pô + cos pp) (sing | z ? x) 2? 
=p (cos ^a cos to $ x eae 2: | sin $ = $ x ) 
+0 (innen; | s H cospsindz, | t Z| | 2? 
-poti ti 
2. The calculations are 
a 9 1 
E = 5 
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Vp — V. (cos pi+sin pj) 


x . y ] 
=V. ! 


| 0 x , 9 y 
ax | yy) oylyi2+y 
u xd xad — x2(x2 + y2)-1/2 | Vx xy? - y2(x2 + y?) 1/2 


x2 + y? x2 + y? 
eee a A 
x? + y? 
1 1 


Vear P 


3. Vx p=0,V-¢=0and 


1 
4. With p constant, the mass of the cone is given by M = pV = 3 7pa°b. By symmetry, we argue that 


R = Zk. Then 
p 
Z= — dV. 
al; 


We can perform the dxdy integration at a fixed z by finding the area of a circular cross-section of the 
. : T . AC az 
cone located at height z above the apex. The radius of this circular cross-section is given by r = + 


Therefore, the volume integral reduces to a one-dimensional integral over z, and we have 


b 2 
-P po 
z=% zn (=) dz 
_ pa? /b* E. 
inpa?b Jo 


_3/P 
pla 
3 
q? 


Solutions to the Problems for Lecture 33 


1. The spherical coordinate unit vectors can be written in terms of the Cartesian unit vectors by 


f =sin0cos pi+sin0sin $ j + cos Ok, 
cos 0 cos ġ i + cos @sing j —sin0 k, 


ô = 
Q — -singi 4 cos $ j. 
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In matrix form, this relationship is written as 


f sin8cos$ sin@sing cos 0 i 
ô | = |cos0cos$  cos0sin$  —sin? j 
ó — sin $ cos $ 0 k 


The columns (and rows) of the transforming matrix Q are observed to be orthonormal so that Q is an 


orthogonal matrix. We have Q-! = QT so that 


i sindcos® cosBcosp —sing f 
j | =|sin@sing cos@sing cos $ ô |; 
k cos 0 — sin 0 0 $ 


or in expanded form 


i=sin0cosp f + cos 8 cos p Ó — sin $ ¢, 
j = sin Osin $ f + cos Osin $ Ó + cos 6 e, 


k = cos 0 f — sin 0 Ô. 


2. We need the relationship between the Cartesian and the spherical coordinates, given by 
x =rsin0cos $, y = rsin@sin®g, z = r cos. 
The Jacobian to compute is 


ox/dr 0Ox/900 0O0x/0$ sindcos® rcos0cos$ —rsin@sing 
dy/dr 0y/00 dy/dd| = |sinÜsin$ rcos@sing | rsin0cosó 


Oz/Oór 0O0z/00 0z/0$ cos 0 —r sin 0 0 
sindcos® cosgcosp —sing 
= r^sinO sinOsin$ cossin p cos $ 
cos 0 — sind 0 


= i? sing (sin? 0 cos? p + cos? 0 cos? p + sin? 0 sin? p + cos? 0 sin? 9) 


— 1? sino (sin? 0 + cos? 9) (sin? $ + cos? p) 


= i? sin 6. 


Therefore, dx dy dz = r? sin dr d0 dọ. 


3. We have 


fav = ud a "ho sin @ drd dọ 
e 
= [ap [ soon [pyar 


= ax [Pa 
0 
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where we have used Int dp =27 and fj sin@ dé = — cos e|» -2 


4. To find the mass, we use the result 


m= fff p(x, y, z) dx dy dz, 
V 


where p is the object's mass density. Here, the density p is given by 


p(r) = po + (pı — go)(r/ R), 


and the total mass of the sphere is given by 
2x pm pR 
mM=/ ll / [po + (pı — po) (r/ R)] r° sin 8 dr do dp 
0 0 JO 


R rà R3 er R3 
= an | por + (Pı — po) 4 dr — 4n E . (1 Lo) 
0 


4 1 3 


The average density of the sphere is its mass divided by its volume, given by 


"ER de 
pP = 700 ras 


Solutions to the Problems for Lecture 34 


1. We begin with 
f = sin 0 cos $ i + sin 0 sin $ j + cos 0 k, 


= cos 89 cos $ i + cos O sin $ j — sin 0 k, 


ô 
Q = -singi + cos $ j. 


Differentiating, 
of n 
a = cos 0 cos $ i + cos0sin $ j — sin 0 k = 6; 
and 
or "NE dcc o E 
36 = —sin0sin $i --sin8cos$j = sind od. 


2. The computations are 


x 129,5 2 Lm 
ASAS Vxf?- 7 
n 1 09,. cos 0 ~ 09 È 
ô = 9) = =? (r) = £; 
rand 3g ein ) r sin0' Yd a r 
er: 69 fcos0 6 
pmo) r Or (r) rsinÜü r` 


A ^ T 
Vagos. VD ind dd 
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3. Using spherical coordinates, for r + O for which 1/r diverges, we have 


1 19 o (1 129 
2 _ 2 = = 
d (=) r2 or (r or C) gs 1130 


Solutions to the Practice quiz: Cylindrical and spherical coordinates 


1. b. We compute using the Laplacian in cylindrical coordinates: 


2/1 19 9 /1 19 1 i3 /1 1 
V = p = p >) = = 3. 
p) pop op \p/ pop p popXp) P 


2. c. When r = xi, the position vector points along the x-axis. Then f also points along the x-axis, Ô 


points along the negative z-axis and $ points along the y-axis. We have (7,6, 9) = (i, — k, j). 


3. c. To find the mass, we use the result 


m= fff p(x, y, z) dx dy dz, 
V 


where p is the object's mass density. Here, with the density p in units of g/cm?, we have 


p =p(r)=10-r. 


The integral is easiest to do in spherical coordinates, and using dx dy dz = r° sin 0 dr d0 d, and com- 


2m pm r5 
m= f n n (10 — r) 2 sin @ dr d0 dp 
0 0 0 


5 
= an | (107? — r°) dr 
0 


E 105 1, 
an ($r y) 


x 3272 g & 3.3 kg. 


puting in grams, we have 


? 31257 
< 3 


Solutions to the Problems for Lecture 35 


1. Let y = y(x). On the one hand, the line integral over a curve C is given by 


[fenis - [reso eras 
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On the other hand, we can parametrize the curve by x — t and y — y(t). Then, 


[rend = [fetu Yan? eva 
C to 

= f ftev» i soa 

= f" fey) eras. 


Solutions to the Problems for Lecture 36 
1. We parameterize a circle of radius R by 
x(0) = Rcos6, y(0) = Rsin6, 
where the angle 0 goes from 0 to 27. The infinitesimal arc length ds is given by 
= „/x' (0 +y(0) de 


= V R2sin2 0 + R2 cos? 0 de 
— R de. 


The circumference of a circle — or perimeter P — is then given by the line integral 


27 
Pee Rd0 = 27R. 
e 0 


Solutions to the Problems for Lecture 37 


1. We integrate u = —yi + xj counterclockwise around the square. We write 


E udr+ | wart f wart | u:dr, 
C C C5 C3 C4 


where the curves Cj represent the four sides of the square. On C; from (0,0) to (L,0), we have y = 0 
and dr — dxi so that / u-dr = 0. On C» from (L,0) to (L,L), we have x = L and dr = dyj so 


Cy 


L 
that n u:dr = n Ldy = L?. On C3 from (L,L) to (0, L), we have y = L and dr = dyi so that 
C2 0 


0 
D u-dr = n —Ldx = L2. The sign of this term is tricky, but notice that the curve is going in 
C3 L 


the -i direction and so is the x-component of the vector field so the dot product should be positive. 
On C, from (0, L) to (0,0), we have x — 0 and dr — dyj so that ri u:dr = 0. Summing the four 


C4 
contributions, we found 
$ u- dr = 212, 
c 


which is twice the area of the square. 
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2. We integrate u = —yi + xj counterclockwise around a unit circle. To parameterize a circle with 
radius R, we write 
x = Rcos6, y = Rsiné. 


Therefore, u = —R sin 9i + Rcos6j and dr = (—Rsin6i + R cos05)d0. We have u - dr = R?d0 and 


270 
pu-dr=[ R? d0 = 2n R?, 
C 0 


which is twice the area of the circle. 


Solutions to the Problems for Lecture 38 


1. We define our coordinate system with the x-axis pointing downward and the origin at the initial 
position of the mass. With 
F = mgt, dr = dxi, 


the work done by gravity as the mass falls a distance h is given by 


h 
w= | F-ar= | mg dx = mgh. 
C 0 


With v; the final velocity of the mass, and with the initial velocity equal to zero, we have from the 


work-energy theorem, 


1 
mgh — ¿mlofl, 


or 


[of| = 4 2gh. 


Solutions to the Practice quiz: Line integrals 


1. d. We have 


ds — V (dx) + (dy)? — TE (dy/dx)? dx = yı + (2x)? dx. 
Therefore, the arc length is given by 


1 
/ v 1-4 4x? dx. 
0 


2. c. We integrate u — —yi + xj counterclockwise around the right triangle. We write 


fuar [ wart | u-dr+/ u:dr, 
c Cy Co C3 


where the curves C; represent the three sides of the triangle. On C4 from (0,0) to (L,0), we have y = 0 


and dr — dxi so that f u dr = 0. On C from (L,0) to (0, L), we parameterize the line segment by 
Cy 
r = (L—s)i+sj as s goes from zero to L so that dr = —dsi + dsj. Therefore, on this line segment, 


u-dr = (x+y)ds = Lds. We have | 


L 
u-dr = n Lds = L?. On Cs from (0, L) to (0,0), we have 
C2 0 
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x — 0 and dr — dyj so that / u: dr = 0. Summing the three contributions, we found 


Ca 
$ ue dr=12, 
E 


which is twice the area of the triangle. 


3. a. Define the x-axis to point vertically upward. The gravitational force is given by F = —mgi and 
the work done by gravity on the way up is -MgXmax and the work done by gravity on the way down 
is M8 Xmax, Where Xmax is the maximum height attained by the mass. The total work done is zero. 


Solutions to the Problems for Lecture 39 


1. 


a) The unrolled cylinder is a rectangle with dimensions as shown on the figure below: 


27a 


The lateral surface area is A = 27tab. 


b) Define the cylinder parametrically as 
r=acos0i+asin0j+zk, forO<z<b and 0<6< 27. 


To find the infinitesimal surface element, we compute the partial derivatives of r: 


or , . . or 
3g asindi+acos07, x^ 
The cross product is 
1 j k 
or or . . : . 
38 xU -asind acos@ 0O| =acosgi+asinó j, 
0 0 1 
so that 
or or 


= Va2cos26 4- à sin? 6 = a. 
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The surface area is given by 


b 2x1 
A= | as= | / adó0 dz = 27tab. 
S 0 JO 


a) The unrolled cone is a circular sector with dimensions as shown on the figure below: 


Va? +b? 


The lateral surface area is the sector of a circle of radius v à? + b2. Its area is found from 


arc length of sector 
circumference of circle 


us x mua? + b?) 


"A ERP 
2 
= na VTT = mab 1 (2) : 


x area of circle 


b) Define the cone parametrically as 


r= " cosbi+ = sind j +zk, for0<z<b and 0<0<21, 


To find the infinitesimal surface element, we compute the partial derivatives of r: 


or az . . az . dr a Pd is 
38 7 sin@i+ 5 cos 0 j, 3z j 0501 „eine Hk. 
The cross product is 
i j k 
or Or . az . aZ . az 
38 x ae —* sing 5 cos 0| = p cos 84 ig, sin 0 j p k, 
¿ cos 0 psinü 1 
so that 
dr or a?z? z 
E 0 20 bam 
30 Oz E cos p sin^ 0 + pi 
2 
_ 92 eee 
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The surface area is given by 


b r2n 2 
az | a 
a=] a=] | —4/1+ 5 d0 dz 
S o Jo b b? 
a a2 b 270 
=-4\ 1+5 d dé 
b +5 f 


= maby/1+ OR 


Solutions to the Problems for Lecture 40 


1. For the paraboloid, we have 
b 
z(x,y) = 5 (2 y!) 


The surface area is given by 


s= f as= | j: | (=) | (E) iav 


Here, 


so that 


$= / [1 + ABS Jat + AP? Jah dx dy 
5 


2b at 
= ah "ES + (x2 + y?) dx dy. 


We integrate in polar coordinates. Let 


x=rcos0, y =rsin0, 


and dx dy = r dr d0. The integral becomes 


2b [?* p^ [gt 
S=7 , n mt” r dr do 


2 4 2 
_ lu /4b 


-7 u’? du 
A a4 /AD2 


_ Amb 2, pore ey 
~ Baz (NS Ir 402 

4 Qi e ay? 
ed ((: | 7) lan) 
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Solutions to the Problems for Lecture 41 


1. Using cylindrical coordinates with the z-axis down the center of the cylinder of length ! and radius 


a, we have for the surface integral on the lateral surface, 
T — ap 4- z£, dS = dS, 
and for the total surface area of the lateral surface, S = 27tal. Therefore, 
EE dS = 2l. 
5 S 
On the top disk of the cylinder, we have 
and for the total surface area of the top disk, 5 — ma”. Therefore, 


l mal 
r-dS = [as= : 
/ 2 Js 2 


The bottom disk of the cylinder will yield the same result, so that 


mal 


pr -dS = 27071 42x = 37a"l, 
S 


which is three times the volume of the cylinder. 


Solutions to the Problems for Lecture 42 


1. To find the mass flux, we use 
r\2 
u(r) = Um (1 - (=) ) k. 


With dS = kr dr d0, the mass flux through a cross section of the pipe is given by 
27 pR r\2 
-dS = 1-—(— dr dO 
[rn [my 
1 
= 2n Rpus, | (1 — s?) sds 
0 
= ¿TER? pl 


This is one-half of what it would be if the entire fluid was moving with velocity um. Using the formula 


nGR* 
fouas- w 


for um, we have 


This result is called the Hagen-Poiseuille equation, which relates the pressure gradient to the radius 


of the pipe for a fixed mass flux. 
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Solutions to the Practice quiz: Surface integrals 


1. d. The parameterization of the torus is given by 
x = (R +rc0s 6) cos $, y = (R4 rcosO)sin $, z — rsin, 
so that the position vector is given by 
r = (R+rcos@)cos~i+ (R 4- rcos0)sin 0j 4- r sinO k. 


The partial derivatives are 


d 
ad = —rsin0cos $i — rsinOsin $j +rcosék, 
0 
35 = —(R+rc0s09) sinpi + (R 4 r cos0) cos pj, 
and 
i j k 
Or or : : : 
3e * T — r sin0 cos $ —r sin O sin $ r cos 0 
—(R+rsin0)sing (R+rc0s0)cos q 0 
= —r(R +r cos 0) cos 80 cos pi —r(R +r cos0)cos8 sinj — r(R +rc0s 0) sin0 k. 
Therefore, 
1/2 
T x 7 = [PR + cos)? cos? 0 (cos? $ + sin? $) + r?(R + r cos 0)? sin? J 
1/2 
= [PR + r cos 0)? (cos? 0 + sin? 2] — r(R 4 r cos), 
and 


dS = r(R + r cos 0) d0 dp. 


2. c. In cylindrical coordinates, u = xi + yj = pp. The cylinder ends have normal vectors 2 and —£, 
which are perpendicular to u. On the side of the cylinder, we have u = Rf and dS = pdS, so that 


pu dS = R [as -zROZRL)-2zR^L 
S 


3. a. We perform the flux integral in spherical coordinates. On the surface of the sphere of radius R, 


we have 
u = zk = (Rcos0)(cos 0f — sin 00), 


and 
dS = f R? sin 8 d0 dp. 
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Therefore, the surface integral over the upper hemisphere becomes 


27 m/2 
[was | n R? cos? O sin 0 d0 dp 
S 0 0 


7/2 1 2 
= an? | cos? 0 sin 0 d0 = zuo | w? dw = de 
0 0 


— R3, 
3 


Solutions to the Problems for Lecture 43 


1. With $(r) = x?y + xy? + z: 


a) Vo = Qxy t y?)i + (x? -2xy)j +k 


b) Using the gradient theorem, fc Vo -dr = $(1,1,1) — $(0,0,0) = 3. 


c) Integrating over the three directed line segments given by (1) (0,0,0) to (1,0,0); (2) (1,0,0) to 
(1,1,0), and; (3) (1,1,0) to (1,1, 1): 


[ veo - [ vorar+ | votre f Vod 
1 1 
=0+ f (1+2y)dy+ [de 
= 3. 


Solutions to the Problems for Lecture 44 


1. With u = (2xy + z2)à + (2yz + x2)3 + (2zx + y?)k: 
a) 
i j k 
Vxu=| 09/0x 9/0y 9/02 
2xy +z? 2yz+x? x+y? 


= (2y — 2y)i + (22 — 2z)j + (2x —2x)k 


b) We need to satisfy 


dp — 2 o — 2 dp — 2 
dy TY Tz, a tx, à 2zx t y. 


Integrate the first equation to get 
$-— few +2?) dx = x^y + x2 + f(y, z). 
Take the derivative with respect to y and satsify the second equation: 


=2yz+x or — = 2yz. 
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Integrate this equation for f to get 
pe f dy = y'z + g(2). 


Take the derivative of p = x?y + xz? + y*z + g(z) with respect to z and satisfy the last gradient 
equation: 
2xz-yw^-s'(z) -2zx--y' or (2) =0. 


Therefore, 2(z) = c where c is a constant, and $ = x*y + y?z + zx +c. 


Solutions to the Problems for Lecture 45 


1. Let Vescape be the magnitude of the escape velocity for a mass launched perpendicular to the Earth's 
surface. When the mass reaches infinity, its velocity should be exactly zero. Conservation of energy 


results in 


1 mM 

5 escape eus R =0, 
or GM 

U sans = MERO — 2gR. 


Therefore, we have 
Vescape — Jf 2g R. 


Solutions to the Practice quiz: Gradient theorem 


1. b. f Vó-dr = $(1,1,1) — $(0,0,0) — 1. 
C 


2. a. Since V x u = V x (yi + xj) = 0, the line integral u around any closed curve is zero. By 
inspection, we can also observe that u = V $, where q = xy. 


3. c. To solve the multiple choice question, we can always take the gradients of the four choices. 
Without the advantage of multiple choice, however, we need to compute $ and we do so here. We 
solve 

op dp 


a =2x+y, a =2y+X, 


2 cd 
ox oy nz 


Oz 


Integrating the first equation with respect to x holding y and z fixed, we find 
p= [exty)dx =x + xy + foa). 


Differentiating $ with respect to y and using the second equation, we obtain 


jt agis Or 3 c2 


dy y 


Another integration results in f(y,z) = y? + g(z). Finally, differentiating $ with respect to z yields 
g'(z) = 1, or g(z) =z +c. The final solution is 


(x,y,z) = x^-Exy--y^ 4- zc. 
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Answer c. is correct with the constant c — 0. 


Solutions to the Problems for Lecture 46 


1. Using spherical coordinates, let u = u,(r,0,p)# + ug(r,0,p)Ó + ug(r,0,#)&. Then the volume 


integral becomes 


2n 1 0 . 1 Qu 2. 
Ji (V -u) dV = ri [ [ E uy) 4 gnomo + e) r^ sin 0 dr d0 d. 


Each term in the integrand can be integrated once. The first term is integrated as 


Fr A (23,7 ru, ) 1? sin6drd6 d$ = Er (s a (ur jar) sinodo 


27 
= / n u,(R,0, 9) R? sin 6 d0 d. 
[U 0 


The second term is integrated as 


ua RE: 2 sin0 dr d0 d aa " 9 (sing d |rdrd 
i LÍ (ein uo) ) sin rasap | PU 59 (sin ug) )rarap 


27 R 
= n kinne en avus Ctr eC: 


since sin (71) = sin (0) = 0. Similarly, the third term is integrated as 


[ET Gara at) son - ara nez p) raras 


-f [ (ug(r, 0,270) — ug(r,0,0)) r dr d0 
= 0, 


since uplr, 0,27) = uplr, 0,0) because q is a periodic variable with same physical location at 0 and 27. 


Therefore, we have 


2m pn 
i (V -u) dV =| f u,(R,0, p)R? sin 6 d0 dp = pu -dS, 
V 0 0 S 


where S is a sphere of radius R located at the origin, with unit normal vector given by 7, and infinites- 


imal surface area given by dS = R sin dé dọ. 


Solutions to the Problems for Lecture 47 


1. With u = x?y ¿+ y2z j +2%x k, we use V-u = 2xy + 2yz + 2zx. We have for the left-hand side of 
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the divergence theorem, 


L pL pL 
I wav=2/ f / (xy + yz + zx) dx dy dz 
V 0 JO JO 


L L L L L L L L L 
=2 / xdx | vay | d+ | ax | yay | zdz+ f xax f ay | zdz 
0 0 0 0 0 0 0 0 0 


= 2(L°/4+ L^ /4 + 15/4) 
EIS. 
For the right-hand side of the divergence theorem, the flux integral only has nonzero contributions 


from the three sides located at x = L, y = Land z = L. The corresponding unit normal vectors are 3, 
j and k, and the corresponding integrals are 


L pL L pL L pL 
puads= f li Lydydz+ | f Dedede + | n I?x dx dy 
S 0 J0 0 J0 0 J0 


—15/24-I9/2-4- I? /2 
= 315/2. 


2. With r = xi + yj + zk, we have V -r = 3. Therefore, from the divergence theorem we have 


frais= [ vorav=3 | av e. 
S V V 


Note that the integral is equal to three times the volume of the box and is independent of the placement 
and orientation of the coordinate system. 


Solutions to the Problems for Lecture 48 


1. With u = f/r, we use spherical coordinates to compute 


1d 1 


u = 2a) = = 


Therefore, for the left-hand side of the divergence theorem we have 


R 
| (v idv =a f DLE 
V 0 


For the right-hand side of the divergence theorem, we have for a sphere of radius R centered at the 


1 AnR? 
huas ETE 7 ArıR. 


origin, dS = f d$ and 


2. With r = xi + yj + zk, we have V - r = 3. Therefore, from the divergence theorem we have 


[ras [vorav=a f av=s (¿nr?) = anr. 
S y V 3 
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Note that the integral is equal to three times the volume of the sphere and is independent of the 
placement and orientation of the coordinate system. 


3. We consider the velocity field given by 


A (xà 4- yj + zk) 


u(x,y,z) — An(x2 + y? + 22)372° 
a) Using spherical coordinates, we have r = rf with r = yx? + y? + 2?. Therefore, 
ENS 
uc amt 


b) We compute V - u for r + 0 using spherical coordinates: 


13 gee ae 
vn zh t-3J) j 


c) We now consider the volume integral of V - u, i.e., 


f yuav. 
V 


If V does not contain the origin, then this volume integral is zero. If V contains the origin, and 


since V - u = 0 everywhere except at the origin, we need only integrate over a small sphere of 


volume V' € V centered at the origin. We therefore have from the divergence theorem 
[Vv -uav = V-udV = u-dS, 
V V! t 
where the surface S' is now the surface of a sphere of radius R, say, centered at the origin. Since 


dS — fdS, we have 
A 
-dS = dS = ^ 
fj, > » Anz R? $, 5 s 


since the surface area of the sphere is A7zt R2. Therefore, 


0, (0,0,0) EV; 
Ir wav = 
V ^, (0,0,0) € V. 


ıYı 


For those of you familar with the Dirac delta function, say from my course Differential Equations 
for Engineers, what we have here is 
V-u=Ad(r), 


where ö(r) is the three-dimensional Dirac delta function satisfying 
ó(r) =0, when r # 0, 


and 
/ ó(r)dV — 1, provided the origin is in V. 
V 


In Cartesian or in spherical coordinates, the three-dimensional Dirac delta function centered at 
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the origin may be written as 


Solutions to the Problems for Lecture 49 


1. The continuity equation as derived in the lecture is given by 


op _ 


Using the vector identity V - (pu) = u- Vo + pV - u, the continuity equation becomes 


Lu Vp+pV-u=0. 


2. We begin with 
d 
5] OA - 3-48 


The divergence theorem applied to the right-hand side results in 


ps-as= | v aav; 
S V 


and combining both sides of the equation and bringing the time derivative inside the integral results 


9p _ 
f (v7) dV =0. 


Since the integral is zero for any volume V, we obtain the electrodynamics continuity equation given 


by 


in 


9p 
3 *V 79 


Solutions to the Practice quiz: Divergence theorem 


1. a. With u = yzi + xzj + xyk, we have V - u = 0. Therefore, 


fuis [iv av =o. 


2. d. 
pras= | (Vnav=3 | av «sni. 
S V V 
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3. d. Computing the divergences, we have 


. 15. 
v: [svi wil =y-y=0 


V-[(1+x)é+(1—y)j] =1-1=0, 
vV. " Mol 24) ] = (2x — y) + (y =2x) =0, 
V-[(x4 xy] =2(x +y) -22-y) = 4y. 


Solutions to the Problems for Lecture 50 


1. With u = —yi + xj, we use ðu2/9x — ðuı/ðy = 2. For a square of side L, we have for the left-hand 


/ (= st) da =2 f dA = 212. 
A ox dy A 


When the square lies in the first quadrant with vertex at the origin, we have for the right-hand side of 


side of Green's theorem 


Green's theorem, 


L 0 0 L 
jon dx + ua dy) = | ode f (-L)ax + f ody + f Ldy =21?. 
c 0 L L 0 


2. With u = —yi + xj, we use dun/dx — du: /dy = 2. For a circle of radius R, we have for the left-hand 


/ E e) dA =a dA = 2nR?. 
A\ Ox dy A 


For a circle of radius R centered at the origin, we change variables to x = R cos0 and y = R sin 6. Then 


side of Green's theorem, 


dx = —Rsin@ and dy = Rcos 0, and we have for the right-hand side of Green's theorem, 


27 
$ a dx + uo dy) = dv dx+xdy) = T (R? sin? 0 + R? cos? 8)d0 = 2n R?. 
C C 0 


Solutions to the Problems for Lecture 51 


1. Let u = ui(x,y,z) i + uo(x, y, z) j + ua(x,y,z) k. Then 
duz dm\. du, ðuz\. du» ðu 
MOOR (S ez dx)? ox oy 


a) For an area lying in the y-z plane bounded by a curve C, the normal vector to the area is i. 


Therefore, Green's theorem is given by 


Qua Qu = i 
[ ( a ) dA fondy + uadz); 
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b) For an area lying in the z-x plane bounded by a curve C, the normal vector to the area is j. 
Therefore, Green's theorem is given by 


du Aug 
dA = d dx); 
LG =| Bu z + udx); 


The correct orientation of the curves are determined by the right-hand rule, using a right-handed 


coordinate system. 


2. We have u = —yi + xj. The right-hand side of Stokes’ theorem was computed in an earlier problem 
on Green’s theorem and we repeat the solution here. For a circle of radius R lying in the x-y plane 
with center at the origin, we change variables to x = R cos $ and y = R sin $. Then dx = —R sin $ and 
dy = Rcos $, and we have for the right-hand side of Stokes’ theorem, 


27 
$ u-dr = P 0n dx + ta dy) = d ydx + x dy) =| (R? sin? p + R? cos? $)d0 = 2x R2. 
C C C 0 


The left-hand side of Stokes’ theorem uses 


i j k 
V xu = |ð/ðəx d/dy 0/0z| = 2k; 
—y x 0 


so that with dS = FR? sin 0 d0 d$, we have 


277 m/2 
x w-as=2r f / k- ê sin 0 d0 dọ. 
S 0 0 


With 
k = cos0 f — sin 0 Ô, 
we have 
k- f = cos 0; 
and 


270 m/2 
[vxo as-o f ap f cos 0 sin 0 d0 
S 0 0 


n/2 


"A 2 
0 = 2nzuR*. 


= 2n R? sin? J 


3. We consider the velocity field given by 


T (—yi+xj 
u = E 
2n \ x2 +y? 
a) Polar coordinates are defined by 


x-—rcosÜ, y=rsin0, i = cos0f — sinb, j = sin 0f + cos 0Ó. 


Substituting into u, we see that 


b) 


c) 
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Let w = V x u. For r #0 and u, = 0, ug = T/27tr, we have using polar coordinates, 


„(19 1 du, (19 /T 
ESSE k (Luo) ras) E (zz) à 


We now consider a surface integral of the vorticity field over a surface in the x-y plane containing 


nm 
S 


Because w = 0 everywhere except at r = 0, we can reduce this integral to be that over a circle of 


the origin, i.e., 


radius R centered at the origin. Then applying Stoke's theorem, 


RN " ?n / TO 
Je kas = |V xu) kas =f u-dr= | (zur) ôR =r. 
S s E o NX2zR 


Since w equals zero every except at r = 0, and its integral over any surface in the x-y plane 
containing the origin is equal to T, we identify w with T times the two-dimensional Dirac delta 
function, i.e., 

w=Td(r). 


This is called a two-dimensional point vortex. The two-dimensional Dirac delta function can also 
be written in Cartesian or polar coordinates as 


1 
&(r) = é(x)8(y) = zêl) 
Solutions to the Practice quiz: Stokes' theorem 
1. b. With u = —yi + xj, we have V x u = 2k and we use Stokes’ theorem to write 


fouearo [( 48-2 f 44 = HR, 
C S 2 


1 
where we have used dS = kdA and the area of the quarter circle is UTR. 


2. 


C, 


With u = 


| 73, one can show by differentiating that V x u = 0 provided 


—y i 
x2 + y? x2+y 


(x,y) + (0,0). However, the integration region contains the origin so the integral is best done by 


applying Stokes’ theorem. We use cylindrical coordinates to write 


— M -—— 
Xy yop 


Then, 


^ 


" l u 27 $ (à u 27 u 
[v rw-as- gu ar= f (2) (dp dp) = A d$ = 2n. 


8. c. With u = —x?yi + xy?j, we have V x u = (x? + y?)k. Therefore, 


1 pl 1 1 1 1 
pu-dr=[(vxu)-as= | i GP ey!) dedy = | za | dy+/ dx | Pas, 
C s o Jo 0 0 0 0 3 
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where we have used dS = k dx dy. 


Solutions to the Problems for Lecture 52 


a) The Navier-Stokes equation and the continuity equation are given by 
9 1 
art (e Vje = cp Nu, V-u=0. 


Taking the divergence of the Navier-Stokes equation and using the continuity equation results in 


V-((u: Vju) = Wy, 


p 
Now, 
d du; \ _ ðu; OU; 
A Ox; («Zt Ox; OX; 
Therefore, 
Qu; du 
2y =p e 
ve Pox Ox; 


b) Taking the curl of the Navier-Stokes equation, and using w = V x wand V x Vp = 0, we obtain 


= V x (u: V)u = vw. 


To simplify the second term, we first prove the identity 


1 
ux(Vxu)= 5; V(u:u) - (u: V)u. 
We prove by considering the ith component of the left-hand side: 
du du 
[ux (V x u)]; = €ijeljfm a = Eiji 
m ag „ 2'i 
Ox] B 1 Ox; J Ox; 


u 19 (u;u;) y i 
20x; 7] Fox, 


2 (9i10jm = Simdj1)Uj 


= viu) - (u V)u];. 


Therefore, using the facts that the curl of a gradient and the divergence of a curl are equal to 
zero, and w = V x wand V - u = 0, we have 


V x(u:Vju-2 Vx (¿Vu 0) -ux (v) 
— —V x (uxw) 
= - [u(V w) — (Vu) + (w: V)u— (u: We] 
= —(w- V)u t (u: V)w. 
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Putting it all together gives us the vorticity equation, given by 


So + (u- V)w = (w- V)u - vV?w. 


Solutions to the Problems for Lecture 53 


1. The electric field from a point charge at the origin should be spherically symmetric. We therefore 
write using spherical coordinates, E(r) — E(r)f. Integrating Gauss's law over a spherical shell of 


radius r, we have 
$5 -dS = E(r) $ dS = Anr?E(r) = L. 
S S £0 


Therefore, the electric field is given by 
E(r) = Lp 


= f. 
Aregr2 


2. The magnetic field from a current carrying infinite wire should have cylindrical symmetry. We 
therefore write using cylindrical coordinates, B(r) = B(p)®. Integrating Ampere's law over a circle 


of radius p in the x-y plane in the counterclockwise direction, we obtain 


$ B-dr = B(o) d dr = 2npB(p) = pol, 
C C 


where I is the current in the wire. Therefore, 


